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Abstract. The rigorous microscopic theory of equilibrium crystal shapes has made 
enormous progress during the last decade. We review here the main results which have 
been obtained, both in two and higher dimensions. In particular, we describe how the 
phenomenological Wulff and Winterbottom constructions can be derived from the mi- 
croscopic description provided by the equilibrium statistical mechanics of lattice gases. 
We focus on the main conceptual issues and describe the central ideas of the existing 
approaches. 
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Figure 1. The free energy of the crystal- vapor interface is given by the integral of the 
anisotropic surface tension Tp over dV . Tt^ d ^ 1 ^ is the (d — l)-dimensional Hausdorff 
measure. 



Part 1. Introduction 

1.1. Phenomenological Wulff construction 
1.1.1. Equilibrium crystal shapes. The phenomenological theory of equilibrated crys- 



tals dates back at least to the beginning of the century | Wu|. Suppose that two different 
thermodynamic phases (say crystal and its vapor) coexist at a certain temperature T. 
Assuming that the whole system is in equilibrium, in particular that the volume v of the 
crystalline phase is well defined, what could be said about the region this phase occupies? 
Of course, the issue cannot be settled in the language of bulk free energies - these do not 
depend neither on the shape, nor even on the prescribed volume v of the crystal. Instead, 
possible phase regions are quantified by the value of the free energy of the crystal-vapor 
interface, or by the total surface tension between the crystal and the vaporQ. Equilibrium 
shapes correspond, in this way, to the regions of minimal interfacial energy. This is an 
isoperimetric-type problem: The surface tension Tp (where, throughout the article, (3 de- 
notes the inverse temperature, (3 = 1/T) is an anisotropic function of the local direction of 
the interface. Thus, assuming that the crystal occupies a region V C M d , the correspond- 
ing contribution Ws (V) to the free energy is equal to the integral of Tp over the boundary 
dV of V (Fig. D. 

The Wulff variational problem could then be formulated as follows: 

(WP)„ Wp(V) — ► min Given: vol(V) = v 

As in the usual isoperimetric case (WP)„ is scale invariant, 

Vs > 0, Wp {d(sV)) = s d - l Wp (dV) . 

Consequently, any dilatation of an optimal solution is itself optimal, and one really talks 
here in terms of optimal shapes. 

The canonical way to produce an optimal shape is given by the following Wulff con- 
struction (Fig. ^): Define 

K,= p| |xGlR d : x ■ n^r^ra)} = f] H p {n) . (1.1.1) 
x In this review, our point of view is that of mathematical physics; for an exposition of the problem from 



the viewpoint of theoretical physics, we refer to RW and references therein 



Figure 2. Function 773(7?) (left) with three half-spaces H(fti), H(n2) and H(n 3 ) (for 
better visibility, only H{n\) has been shaded). The intersection of all such half-spaces 
gives rise to the corresponding Wulff shape (right). 



It would be convenient to normalize K, as 

V vol (^) 

We refer to tC\ as to the normalized, or unit volume, Wulff shape. The variational theory 
of (WP)^, which we briefly address in the subsequent subsection, states that any solution 

to (WP) t can be obtained by a shift of the corresponding dilatation K v = ^vK,\ of K,\. 

1.1.2. Variational methods. The corresponding literature is rather rich and diverse, 
here we merely attempt to facilitate the orientation of the reader and to introduce some 
notations which will be useful in the sequel. 

Since the half-spaces Hp (n) in (1.1.1) are convex, so is the Wulff shape K.. Furthermore, 
in all the problems we consider here, the surface tension 773 is bounded above and below, 

< min Tfl(n) ^ max r«(n) < 00. (1-1.2) 

Accordingly, equilibrium crystal shapes are bounded and have non-empty interiors, G 
int(/C„) . 

The fact that K, is optimal follows from the general Brunn-Minkowski theory: Let r^* 
be the support function of JC, Tg*(x) = sup{y • x \ y G /C}. Of course, if the homogeneous 
extension of 773 

A 1 X \ 

Tp{x) = \\x\\ 2 rp ( "pFjpJ ' (1.1.3) 

is convex, then Tp and t%* coincide. In general r^* is the convex lower-semicontinuous 
regularization of r^, in particular rj* ^ 773. Nevertheless, for the Wulff shape /C, 

W?(K) = I T*p*(n x )dHi d -V = [ T^ftJdHj?- 1 ). 

JdK JdK. 

where, as before, n x is the outward normal to dV in x and 7-^ d_1 ) is the (d— 1) dimensional 
Hausdorff measure in R d . 
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On the other hand, the action of the regularized functional WS* could be extended to 
any compact set V C M. d in terms of the mixed volume 

W*p*{V) = liminf-(vol(V + e£) -vol(V)), 

the latter definition coincides with the integral definition of W^* for regular V . The 
Brunn-Minkowski inequality [}5ch| 



d 



vol{A + B)^z [vol(A)d +vol(S)3 

implies that for any regular V with vol (V) = vol (/C), 

Wp (V) > (V) ^ dvolQC) = WpQC). 

Of course, we have been rather sloppy above, and we refer the reader to the works 
[Ta|, |Fj and [FM for the comprehensive discussion and results, including the history of 



the variational Wulff problem. The language employed in the latter works is that of the 
geometric measure theory, and we proceed with setting up some of the corresponding 
notation which will also turn out to be useful for the Li-approach to the microscopic 
justification of the Wulff construction, as described in Part 2 of this review. In the latter 
case, the macroscopic state of the system will be determined by the value of an order 
parameter which specifies the phase of the system. In the systems that we will consider, 
the pure phases are characterized by their averaged density, which are encoded by two 
values pi(/3) and Ph(P), for example ph for the crystal and p\ for the vapor. (In fact, 
we shall derive all the results in the symmetrized spin language, in which case the two 
values will be ±m*(/3), where m*{0) is the spontaneous magnetization (see Section 2) at 
the inverse sub-critical temperature > f3 c ). For a given temperature, it is convenient 
to replace this order parameter by a parameter with values ±1. We suppose that the 
macroscopic region of M. d where the system is confined is the unit torus T d = (M/Z) rf . 
The macroscopic system is described by a function v taking values ±1 and the fact that 
v r = 1 for some r in T d means that locally at r the system is in equilibrium in the phase m* . 

For any measurable set V in T d , the perimeter of V is defined by 

V(V) = su p|y div4>(x)dx | <f) G C 1 (T d ,W d ), \<f>\ < lj . (1.1.4) 

A function v with values ±1 is said to be of bounded variation in T d if the perimeter of the 
set {v = 1} is finite. We denote by BV(T d , {±1}) the set of functions of bounded variation 



in T with values ±1 (see [ EG | for a review). For any v in BV(T , {±1}), there exists a 



generalized notion of the boundary of {v = 1} called reduced boundary and denoted by 
d*v. If {v = 1} is a regular set, d*v coincides with the usual boundary dv. Furthermore, a 



blow-up Theorem (see [ EG | p. 199) ensures that for all x in d*v an approximate tangent 
plane can be defined locally. This will imply the existence of a unit vector n x called the 
measure theoretic unit normal to {v = 1} at x. For any x in M d and any vector ft, we 
define the half spaces 

H + (x,n) = {y G R d | (y - x) ■ n > 0} , 
H~(x,n) = {y G R d \ (y - x) ■ ft < 0} . 
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H (x,fi) 



d*v 



{« = !} 



Figure 3. Measure theoretic unit normal to {v = 1} at x 



Then for all x in d*v , there is a unit vector n x such that 
lim -^vol (B(x, r) f]{v = l}f] H + {x, 

& ^ Vd (S(ar,r)n{« = -l}n ff ~( a? > 



I?. 



0. 
0. 



where B(x,r) is the ball of radius r centered in x. The previous property shows that the 
reduced boundary is not too wild (see Fig. |3|). In fact, it is possible to prove that a set of 
finite perimeter has "measure theoretically a C 1 boundary" . 



The functional Wp can be extended on Li(T rf 



-ttI) as follows 



Wp{v) 



if u e BV(T d ,{±l}) 
otherwise. 



(1.1.5) 



Under the assumption that the homogeneous extension ( |L1.3|) of rp is convex, a result by 
Ambrosio and Braides (see [ AmBi], Theorem 2.1) ensures that Wp is lower semi-continuous 
with respect to Li convergence. In certain cases (attractive interactions) the convexity of 
Tp can be derived from the properties of the corresponding microscopic system as will be 
explained later. 

To any measurable subset A of T d , we associate the function 1a = ^-A c — 1a and simply 
write Wp(A) = W/j(]Ll). In this new setting, the isoperimetric problem is to find the 
minimizers of 



mm{Wp(v) | v G BV(f d ,{±l}), 



m* v r dr\ < m}, 



(1.1.6) 



where m belongs to ]m({3) , m* ((3)[. The parameter fh is chosen such that the minima of 
the variational problem above are translates of the set fc m deduced from the Wulff shape 
K, by dilatation in order to satisfy the volume constraint. This restriction enables us to 
exclude pathological minimizers which occur from the periodicity. Nevertheless, notice 
that the precise shape or the uniqueness of the minimizers of the variational problem will 
be irrelevant for the microscopic derivation of the Wulff construction. 

1.1.3. Stability properties. In two dimensions Wulff solutions to (WP) are stable in 
the metric of Hausdorff distance: let V be a connected and simply connected subset of IR 2 
with a rectifiable boundary dV. Assume that Area(V) ^ 1. Then, 



mmd H (V> + /Ci) < c^WpW-WpilCi 



(1.1.7) 
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This result has been established in [DKS] as a generalization of the classical Bonnesen 
inequality. 

If V consists of several connected and simply connected components, V = V™ =1 Vi, and 
the total surface tension of V is close to the optimal, 

n 

i=i 

then, again assuming that Area(V) = ^27=1 Area(V^) 1, an easy consequence of ( |1 . 1 . 7[ ) 
implies (see (2.9.7) and (2.9.8) in [DKS]) that actually all but one components of V are 



small, and that the only large component, say V±, is close to a shift of K\. Namely 

n n 

^Area(^) < c 2 e 2 and J^W/j^) ^ C3£ > 



i=2 i=2 



and V% satisfies ( |1.1.7| ). 

These stability properties are indispensable for a sharp justification of the phenomeno- 
logical Wulff construction directly from the microscopic assumptions on the local inter- 
particle interactions (see Section ^5| of Part ||) . 

As far as we understand, stability properties of higher dimensional isoperimetric prob- 
lems are much less studied. Already in three dimensions the Hausdorff distance is, of 
course, not an adequate measure of stability. Trivial rate-free stability properties in Li 
simply follow from the uniqueness of Wulff solutions and the compactness of BV-balls 
in Li. On a more qualitative side there are well studied stability properties in the class 



of convex sets fSchfl and, also, for sets with a smooth boundary [Ha|. We feel, however, 
that the statistical stability under the microscopic approximations in the problems we 
consider here might be better than the impartial stability of the corresponding variational 



problems. A result of this sort is supposed to appear in [BIV|. 



1.1.4. Winterbottom problem. The Wulff variational problem provides a description 
of an equilibrium crystal shape deep inside a region filled with gas phase. If, however, the 
spatial extent of the system is finite, it may happen that the boundary of the surrounding 
vessel exhibits a preference toward the crystal phase. In such a situation, the equilibrium 
state may not be given by the Wulff shape anymore, but may have the crystal attached 
to the boundary. We discuss briefly the simplest model of such an interaction between an 
equilibrium crystal and an attractive substrate. Suppose, for simplicity, that our system 
is contained in the half-space H = {x G R rf : x(d) 0}; the boundary of this half-space, 
the hyperplane It) = {x G M. d : x(d) = 0} represents the boundary of the vessel and is 
called the wall. We also suppose to simplify the analysis, and because these assumptions 
will always be satisfied, that rp{n) = 173 (— n), and that the homogeneous extension of rp 
is convex^]. 

To model the degree of attractiveness of the wall, we introduce a new thermodynamical 
quantity, the wall free energy r bd (/3, 77), which depends on both the inverse temperature j3 
and the "chemical structure" of the wall 7], and modify the free energy functional accord- 
ingly, 

Wp,v(V) = Wp(Y) + (r bd (A v) - t$ (dv n ») , 



2 In the models we consider in this paper, this is a consequence of FKG inequality. 
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Figure 4. The Winterbottom shape is obtained by taking the intersection between the 
Wulff shape and the half-space {x(d) ^ — Tbd(/3, ??)}, and rescaling the obtained body. 



where = Tp(ed), G M. d with ed(&) = ^d- The wall free energy replaces therefore the 
surface tension Tp along the wall. At equilibrium, a thermodynamical stability argument 
shows that r bd (/3, rj) ^ (this can also be proved in some microscopic models, see Part |j), 
so that this last term is always non-positive. The new variational problem is 

(WBP)„ W/9,„(V) — ► min Given: V C H, vol(V) = v . 

It has first been studied in [Wi] and is called the Winterbottom variational problem. Let 
us now discuss what its solution looks like. It turns out that there are three cases to 
consider: 

1- T bd (A?7) =T* p 

In this case, Wp^iV) = Wp(V) and therefore the solution is the Wulff shape as- 
sociated to Tp. The equilibrium crystal is not attached to the wall. This can happen 
even if a priori the chemical structure of the wall is such that it is energetically favor- 
able for the crystal to lay on the wall, see Part |3] for a discussion from a microscopic 
point of view. 

2. \r bd {M\ <r*p 

Now the wall is really attractive for the crystal shape. The solution of the varia- 
tional problem is given by a suitably rescaled version of the following set (see Fig. |j) , 

IC W = K n {x g R d : x(d) ^ - r bd (/3, 77)} 

so that the volume constraint is satisfied (notice that this variational problem is still 
scale invariant); see [KP] for a simple proof. 

3. r bd (/3,??) = -Tp 

This is a somewhat pathological case. Indeed, the solution of the variational 
problem is completely degenerate, the solution being unbounded. A minimizing 
sequence is, for example, 

R n = {x G H : \x(k)\ ^ n, k = 1, . . . , d - 1, x(d) < n l ~ d v} . 

As n — » 00, R n covers the whole wall with a film of vanishingly small width; the 
limiting value of the surface free energy functional is 0. This describes the regime of 
so-called complete wetting where the wall so strongly prefers the crystal that it wants 
to prevent any contact with the gas phase. 
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1.1.5. Microscopic justification. Microscopic models we consider here are simple lat- 
tice gas type models (in the magnetic interpretation), which are going to be defined pre- 
cisely in the next section. The prototype situation when the Wulff construction is thought 
to be recovered as a law of large numbers as the size of the microscopic system tends 
to infinity could be loosely described as follows: Suppose that the particles of a certain 
substance live on the vertices of the integer lattice so that each vertex of Z d could be 
either occupied by a particle or remain vacant. Thus, various particle configurations n 
could be labeled by points of {0, l} Ijd , where one puts n$ = 1 if there is a particle at site 
i G Z d , and rtj = 0, otherwise. These random configurations are sampled from a Gibbs 
distribution P, which takes into account the assumptions on the microscopic interactions 
between the particles. The strength of the interaction is quantified by the value (3 = l/T 
of the inverse temperature; the larger (3 (respectively the smaller the temperature T) is, 
the stronger is the interaction. In many instances sufficiently low temperatures give rise 
to two stable phases - the low density phase (which we call vapor) with an average parti- 
cle density per site pi and the high density phase (crystal) with a corresponding average 
density p h , < pi < p h < 1. 

Suppose now that all the particles are confined to a large finite volume vessel Ajy C 
Z d , where the subindex N indicates the linear size of A^; we put for simplicity |Ajv| = 
N d . Let us fix p G (pi,Ph) an d ask what are the typical geometric properties of particle 
configurations n under the conditional measure P ^ ■ | ^ieAjv Jij = pN 1 ^ . In other words, 

we fix the total number of particles pN d in such a way that it falls in-between the two 
stable values piN d and PhN d . 

The prototype law of large numbers result we have in mind is schematically: 




\ 



i£A N 



1. 



Thus, with an overwhelming P ^ • | ^2 ieAN rii = -probability particle configura- 

tions n on An, n G {0, l} Ajv , obey the following phase segregation pattern: A^v splits 



into two regions, A^r = A^ V A l N , where A^ is occupied by the high density phase, and, 
respectively, A^ by the low density one. The relative volume of A^- can be recovered from 
the canonical constraint 



Ph\^%\ 



+ Pi|A{J 



pN a 



and the shape of A^ is asymptotically Wulff. 

There is a long way even towards making the above statement precise - we should define 
the microscopic models, quantify the notion of phases, in particular of phases over finite 
volumes, and explain how the surface tension is produced in the large limit. 
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1.2. Microscopic Models 

1.2.1. Models with finite-range ferromagnetic 2-body interactions. We want to 
introduce mathematically precise realizations of the models discussed in subsection |1.1.5| . 
As described there, our interest lies in models of lattice gases. For simplicity we restrict 
our attention to a particular subclass of such models, which enjoy several nice properties, 
the Ising models with finite-range ferromagnetic 2-body interactions. 

We consider a family of random variables rii, i 6 Z d , taking values and 1. Any site i 
of the lattice Z d is either occupied by a particle, in which case rii = 1, or empty, in which 
case Hi = 0. The random variables rii are called occupation numbers and they completely 
describe a configuration of the lattice gas. We consider a formal Hamiltonian of the form 

\ ^2 Kij n i n j ) 

the 2-body interactions are such that = ifiujik, ^ and Kij = if ||i — > 
r, where r is the range of the interaction. We introduce two parameters, the chemical 
potential fi and the inverse temperature 0, and set A <s Z d . The Gibbs measure in A with 
boundary condition n E {0, 1} Z is the probability measure on ({0, 1} Z ,-4), with A the 
usual product cr-field, defined by 



1 exp(/3// E igA rii + (3 S2 nirij) if = n;, for all z A, 



{ij}nA^0 

otherwise, 



where 

EiGA n i =N 

Two types of boundary conditions are particularly relevant for us, the 1 b.c, corresponding 
to setting n = 1, and the b.c, n = 0. We also need a different kind of boundary 
conditions: The Gibbs measure in A with free boundary conditions is the probability 
measure on ({0, 1} A ,JT A ) defined by 

u A,fi( n ) = —fT~ exp(/3/i ^2 rii + ^ mnj) . 
A,/i ieA {«,j}cA 

These measures describe the lattice gas in the Grand Canonical Ensemble, in which the 
total number of particles, or equivalently the density p{n) = m ^ ieA n«, is not fixed. The 
description of a gas in the Canonical Ensemble corresponds to the conditioned measure 

u A,»,n( ■ I P( n ) = P) > 

with p G Range(p) (this measure is obviously independent of /i). The existence of the Gibbs 
states is^- = lim A y Z d for n = 0, 1 or free, can be easily proved using correlations 

inequalities; moreover, it is unique if fi / — 5 Ej Jbj- Restricting the chemical potential 
to the particular line \i = — | £^ Joj , it can be proved that there exists a critical value 
00 > C > such that 

• For all < C , there is a unique Gibbs state and -J^p) = 1/2. 



For all > C , p ft 09) = ^(p) > 1/2 > ^ (p) = 
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It is rather convenient to work with another, equivalent, formulation of these models, 
in which the symmetries present when ji = — \ Ylj Joj are more transparent; this is the 
magnetic interpretation. To do this, we introduce a new family of random variables Uj, 
i£Z d , defined by 

(jj = 2rii - 1 . 

The random variables Oi therefore take values in { — 1, 1}; o-; L is called the spin at the site 
i. Expressed in these variables, the model is defined through the following Gibbs measure 
in A with boundary conditions a G { — 1, l} zd , 

{-^ — exp(/? ^2 hi <?i + (3 ^2 Jij <7i<Tj) if <Ti = &i, for all % A, 
ieA {iJ}nA^ 
otherwise, 



where hi G R are called the magnetic fields and the coupling constants Jij = J^_j^ satisfy 
Jij ^ and Jij = if \\i — > r. A configuration a such that <jj = <7j, for all i A, 
is said to be compatible with b.c. a in A; the set of all such configurations is denoted by 
^a,ct- We are particularly interested in the + and — b.c. corresponding respectively to 
a = 1 and a = — 1. The Gibbs measure in A with free b.c. is the probability measure on 
({-1,1} A ,^ A ) defined by 

l i A,h( (J ) = -Jj- ex P{P ^2hiai + f3 J ij Wj) ■ 

A,h ieA {«,j}cA 

Expected value w.r.t. these measures are denoted with brackets notations, (')jiah' •■■ 

In the magnetic formulation, the Canonical Ensemble corresponds to fixing the value of 
the magnetization (density) m(a) = jx\^2 i( z\0-i, 

A,*,h( • i m ( a ) = ™) . 

where rh £ Range (m). If hi = h for all i, then the (infinite- volume) Gibbs states fj,^ . for 
+, — and free b.c. can be shown to exist; it is always unique when /i / 0. The phase 
transition statement takes now the following (simpler) form: There exists oo > (3 C > 
such that 

• For all < f3 c , the Gibbs state is unique and (m)^ = 0. 

• For all > ft, m*(f3) = (m)J > > (m)t fi = -m*{0). 

We will use the terminology Ising models to refer to the lattice gases in the magnetic 
formulation. When h = 0, we will generally omit it from the notations. 

Ferromagnetic models are particularly well-suited for non-perturbative analyses. In- 
deed, they enjoy several very useful qualitative properties, most of which taking form 
of correlation inequalities. Of particular importance for us are the following statements 

(°"A = UieA^)'- 

WA<TB)i >h > WA)i th {<TB)i th , 
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provided hi ^ for all i (1st and 2nd Griffiths', or GKS, inequalities |pq , KS|); also, 



& 



2 



for all i, j and £;, provided hi ^ for all / (GHS inequalities | GHS| ); finally 



for any increasingP] functions / and g, and any h E M (FKG inequality [FKG]). Observe 



that any b.c. can be obtained starting with free b.c. and applying suitable magnetic fields 
on the spins on the inner boundary of A, where the inner boundary of a set A C Z rf is 
defined as 

d in A = {ie A : 3j £ A, t ~ j} , 
where i ~ j means that Jj,- 7^ 0. Similarly, we define the (exterior) boundary of A by 

dA = {i # A : 3j E A, j ~ j} . 

1.2.2. 2D nearest-neighbors ferromagnetic Ising model. A particularly simple mem- 
ber of the above-mentioned class of models is the two-dimensional nearest-neighbors Ising 
model, in which Jjy = if i and j are not nearest-neighbors, and Jjj = 1 if they are. This 
model has still additional remarkable features. First, even though this only plays a very 
marginal role in this review, it is the only one for which it is possible to compute explicitly 
various quantities (free energy, surface tension, correlations, ...). Of more importance for 
our purposes is the property of self-duality that it enjoys. 

The nearest-neighbors model admit a geometric description in terms of very simple 
objects, the contours. To define contours in the present context, it is useful to introduce 
the notion of the dual of the lattice Z 2 . The dual lattice is the set of dual sites 

Z 2 = {xEM 2 : x + EZ 2 }. 

To each edge e = (x, y), x, y E Z 2 , we associate a dual edge e* connecting nearest-neighbors 
dual sites, which is the unique such edge intersecting e (as subset of M 2 ). 

Now, if we consider the Ising model in A d Z 2 with b.c. a, a configuration a E Qa,w is 
entirely determined by giving the following set of dual edges, 

{e* : e* dual to e = {i,j), {i,j} nA/0, aiUj = — 1} . 

The maximal connected components of these dual edges, seen as closed line segments in 
M 2 , are called contours. We denote by 7(0") the contours of the configuration a. The 
boundary dj of a contour 7 is the set of all dual sites belonging to an odd number of the 
dual edges composing 7. A contour is said to be closed if $7 = 0, otherwise it is open. 

A set A is Z 2 is simply connected if Ui£a{ x ^ ^ 2 ■ \\ x ~ ^ V^} is a simply 
connected subset of M 2 . 

Given A C Z 2 , its dual is A* = {i E Z 2 : 3j E A, \\j - = 1/2}. A family of 
contours is said to be A* -compatible if they are disjoint (as sets of bonds and sites) and 

3 A function / : { — 1, 1} Z — > R is increasing if f(a) ^ /(c') as soon as <Ji ^ a[, for all i; it is called 
decreasing if — / is increasing. 

4 The fact that this model is sei/-dual is very convenient, but is not required anywhere. What we need 
is to be able to control precisely the dual of the model; for example, the Ising model on the hexagonal 
lattice is not self-dual, but it would be possible to prove the same kind of statements for this model as for 
the one on the square lattice. 
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are included in A*. A family of contours 7 is said to be (A, a) -compatible if there exists a 
configuration a G Oa,^ such that 7(0") = 7. It is easy to show that for simply connected 
A, A*-compatibility of a family of closed contours is equivalent to (A, +)-compatibility. 
The measure Mact can ' 3e easu Y written in terms of these objects; for any a G £1a,ctj 



^A,» = ^7exp{-2/3 Yl HI' "• 2 - 1) 



where | -y | is the number of edges in 7 and 

Z§(A) = £ exp{-2/3^| 7 |}^ £ II^^)- (1-2-2) 

7 (A,CT)-comp. 7G7 7 (A,CT)-comp. 7ST 

We now discuss the property of self-duality. Let A <s 1? be simply connected. We 
consider the model at inverse temperature f3* in the box A* d Z 2 , with free boundary con- 
ditions. There exists another graphical representation for this model, the high-temperature 
representation, which results from writing 

tF"** = cosh/3*(l + OiOj tanh/T) , 

opening all the brackets and expanding. After a simple summation over cr, this yields 

Z^=C(A) (tanh^)S 767 W= C (A) £ ]\w*(r,P*) 

7 A*-comp. 7 A*-comp. 7G7 

= C(A)Z /3 *(A*), (1.2.3) 

where C(A) is some constant which only depends on the set A. Setting tanh/3* = e~ 2/3 , 
we see from flIT2D and ( |L2^ ) that zJ(A) = Z^*(A*), since A is simply connected. In the 



same way, we can expand the 2-point function, for example, and get the following very 
useful identity 

fa"i>£+= E ^(A), (1-2-4) 
where the sum is over all open contours A such that d\ = {i,j}, and 

zp\a*\\)= Yl U w *^n- 

7 closed 7G7 
(7, A) A*-comp. 



Identity (1.2.4) is the so-called random-line representation for the 2-point function of the 



Ising model, and plays a basic role in the approach to the DKS theory of Part H (see 



[ PV2j , [PV3] for much more details on this topic). What is particularly useful is that the 
weights , which we have defined for an open contour, can be immediately extended to 
any family of A*-compatible contours (closed or open). In particular, if 7 is a family of 
A*-compatible closed contours, then the following identity holds 

^(7)=< + (7C 7 (-))- 

Applications and further results about the random-line representation are given in Sec- 
tion |3.4| and in Part ||. The results stated above also hold when the coupling constants 
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are allowed to vary from edge to edge, provided they remain ferromagnetic; if we denote 
by J(e) the coupling constant at edge e, then the duality relation takes the form 

tanh(/3*J*(e*)) = e" 2/3J(e) . (1.2.5) 

1.2.3. Kac models. In the original van der Waals Theory, the occurrence of phase transi- 
tions is due to long range attractive forces between molecules. In its statistical mechanics 
formulation, these forces are described by Kac potentials that depend on a positive scaling 
parameter e which controls the strength and the range of the potential (see [ KUH|1 ). The 



first probabilistic approach to this model was made in the celebrated paper of Lebowitz 



and Penrose | LePe j. 



In dimension d, Ising systems with Kac potentials are defined by Gibbs measures with 
potentials depending on a scaling parameter e > 

Vi,jeZ d , J? tj = e d J(s\\i-j\\ 2 ) , 

and J is a non-negative, smooth function supported by [0, 1] and normalized so that 

dr J(\\r\\ 2 ) = 1. 

i 

The Gibbs measure on the domain A is denoted by \x e A . The constant e will be so that 
the system has finite but long range interaction. It is convenient to consider interaction 
parameters of the form e = 2~ m (m is typically assumed to be large but fixed). 

This model bridges the finite range models and the mean field models. In particular, 
if the range of the interaction, i.e. e _1 , is scaled proportionally to the number of spins 
then the statistical properties of the system can be recovered from a mean field functional. 
In the true thermodynamic limit, when e is kept fixed while the number of spins goes to 
infinity, the behavior of the system cannot be described by the mean field continuum limit. 
Nevertheless, by localizing in finite size regions it is possible to derive some informations 
from the mean field functional. This strategy was used to recover the phase diagram of 
the model and to prove that it is arbitrarily close to the one of the mean field model when 
e goes to 0. More precisely, let us recall the following result which has been proven by 



Cassandro, Presutti [CP and by Bovier, Zaharadnik |B2fl (see also 



Theorem 1.2.1. For any (3 > 1, there is £q > such that for any e smaller than Eo a 
phase transition occurs and there are at least 2 distinct pure phases [ij and \x~ . 

If (5 > 1 , there is a breaking of symmetry and the spontaneous magnetization is denoted by 
fj,f(ao) = m*. Define m* = \im E ^om*. This Theorem was proven via a renormalization 



procedure which we shall describe in Subsection 2.3.1. 



1.2.4. Surface tension. We fix n a vector in and consider an orthonormal basis 
(e*i, . . . , &d-i, n). Let A(iV, M) be the parallelepiped of M. d centered at with side length 
N for the sides parallel to (e±, . . . , ed-i) and side length M for the sides parallel to n. The 
microscopic counterpart of A(N,M) is denoted by A(iV, M). The boundary dA(N,M) is 
split into 2 sets 

dtA(N, M) = {i£ dA(N, M) \ i.n ^ 0}, 
dZA(N, M) = {i£ dA(N, M) \ i.n < 0}. 
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+ + + + 




Figure 5. Definition of the surface tension. 



We fix the boundary conditions outside A(N, M) to be equal to 1 on d^A(N, M) and 

to —1 on dZA(N,M). The corresponding partition function on A(N,M) is denoted by 
7 p 

Zj A(Af,M),n,±- 

Notice that any configuration a contributing to the partition function Z^^y M ^ - ± con- 
tains a ±-contour 7 which crosses A(N, M) under the "averaged" direction orthogonal to 
n (Fig. |). Such a contour is absent in the configurations a contributing to partition func- 
tions Z^/^y M % + with pure boundary conditions on dA(N, M). This contour represents the 
microscopic ±-interface under the direction H. 

Definition : The surface tension in the direction H £ S^ -1 is defined^ by 

rp(n)=]ixn lim - * log ^f^± . (1 . 2 . 6) 

□ 

The proof the existence of the surface tension can be found in many papers ( [ |Ab2 1 , [Pf2] 



to mention a few). A general approach has been developed by Messager, Miracle-Sole and 
Ruiz fMMRll . The core of their proof is the sub-additivity of the sequence of finite- volume 



approximation to (ft) which is obtained by means of FKG inequality. The proof is also 
valid for a wide range of models like Ising models with finite range interactions, Potts and 
SOS models. Furthermore, they showed that surface tension can be defined with paral- 
lelepipeds A(N, M^r), where Mn is a function of N which diverges as N goes to infinity. 
More general domains can also be considered provided they contain a parallelepiped of 
the type A(N, M N ). 

The convexity of the homogeneous extension of rg (see ( |1.1.3|) ) is a consequence of the 
pyramidal inequality proven in Theorem 3 of RMMRfl : Let A), • • • , A d be d + 1 points of 



l d and denote by (Aj)j<^ the simplex defined by these points. Let Hi be the unit normal 



^Notice that surface tension is sometimes defined with an extra multiplicative factor 
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to Aj and |Aj| its area. Then, the pyramidal inequality says 

d 

|A | Tfs(n ) ^ ^2 \ A i\ r /3(^)- 
1=1 

Note also that the homogeneous extension of tr is continuous because it is locally 
bounded and convex. Furthermore, tr is uniformly positive on S rf_1 . This follows from 
the fact that the surface tension r^(no) in the direction Hq = (1, 0, . . . , 0) is strictly positive 
as (3 is larger than f3 c (see Lebowitz and Pfister [ LePf| ]). 



1.3. Scope of the theory 

The key notion behind the attempts to give a rigorous meaning to the type of the phase 
segregation phenomena, which have been vaguely discussed in Subsection |1.1.5| , is that 
of renormalization or coarse graining. The energy (probability) competes with the 
entropy (number) of microscopic configuration in the corresponding energy shells. Macro- 
scopic quantities like surface tension are produced in the aftermath of the entropy/energy 
cancelation, which is to say that in order to derive large- iV (iV-linear size of the system) 
asymptotics one should renormalize appropriate microscopic objects. The appropriate ob- 
jects here are, of course, microscopic phase boundaries, which decouple between different 
"large" microscopic phase regions. These renormalization procedures could follow two dif- 
ferent trends, depending on whether the renormalized (mesoscopic) structures keep track 
of the microscopic or macroscopic state of the system. 

1.3.1. Dobrushin-Kotecky-Shlosman Theory. The coarse graining of the DKS theory 
closely follows microscopic phase segregation patterns. Basic tools comprise a fluctuation 
analysis of the microscopic phase boundaries and sharp uniform local limit estimates 
over domains encircled by such boundaries. Thus, the notion of finite volume phases 
is quantified by the rate of the relaxation of the statistics of microscopic observables inside 
the microscopic phase regions towards the corresponding equilibrium values. 

The theory has been developed using the low-temperature cluster expansions in the 
seminal monograph [DKS]. Our exposition in Part 3 is non-perturbative and follows the 



works |Pf|, 0, Pi, TPV| , [gcS| and |S§. By and large the existing results are confined 



to the simplest two-dimensional models (percolation and nearest neighbor Ising). 

1.3.2. Li-Theory. The renormalization approach of the Li-theory is, in a sense, opposite 
to that of DKS. In the latter case the principal coarse grained objects (skeletons, see Part ||) 
are built upon underlying families of large microscopic contours. Such information is 
waved out in the Li-approach, and the basic renormalization objects here are the local 
(mesoscopic) order parameters or, in the spin language, locally averaged magnetization 
on various length scales. The idea is that on sufficiently large scales local averages of the 
magnetization are, with an overwhelming probability, close to one of the two equilibrium 
values ±m*. Thus, under the renormalization, configurations are characterized by their 
phase labels on different mesoscopic blocks. The objective of the Li-theory is to describe 
typical mesoscopic magnetization profiles (or their phase labels) under a relaxed canonical 
constraint of shell type. Unlike in the DKS case, the mesoscopic phase labels are classi- 
fied by their proximity to various macroscopic states. Combinatorial complexity of this 
approximation is reduced by an exponential tightness property of the mesoscopic phase 



labels (for a general claim of this sort see Theorem 2.2.1), which enables to restrict atten- 



tion only to Li-compact subsets of feasible macroscopic states, namely to the phase-sets 
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of finite perimeter. The core of the compactness estimates is based on the renormaliza- 
tion decoupling techniques introduced in |Pil[| an d on the methods developed to control 
the phase of small contours by [O], [PV2], [pcS3|] and | ISc ] . These techniques are robust 



enough to be applied on a renormalized scale in any dimensions in a non perturbative 
setting. 

Our exposition in this review is based on the work of [Bl| with, though, one exception 
- we specifically stress that all the relevant estimates of the Li-theory are obtained on 
appropriate finite scales. The validity of Lemma 2.4.1 up to the slab percolation threshold 
follows from the results of [CePi 



1.3.3. Boundary Phenomena. Parts || and y provide a derivation of Wulff construction 
from the basic principles of Equilibrium Statistical Mechanics. Part Qis concerned with a 
study of the effect of the boundary conditions on the macroscopic geometry of the phase 
separation. In particular, it is shown how the interaction with the boundary of the vessel 
can be analyzed, and used to provide a derivation of Winterbottom construction. The 
relationship between the macroscopic geometry in this case and the wetting transition is 
also discussed. The presentation follows [PV2] for the 2D case, and [BIV] for the higher- 
dimensional ones. 



1.3.4. Bibliographical review. The rigorous investigation of the macroscopic geometry 
of phase separation under a canonical constraint certainly started with two seminal pa- 
pers of Minlos and Sinai in 1967-68 [MSI, MS2 |. In these papers, the authors considered 
nearest-neighbor very low temperature Ising models in arbitrary dimensions d ^ 2, even 
though they only wrote down the proof explicitly in the case d = 2. Their results could 
be roughly stated in the following way: At sufficiently low temperatures, typical configu- 
rations of the Ising model in the exact canonical ensemble over finite vessels of linear size 
N, consist of a single large contour whose shape is "nearly a square" , whereas the rest of 
the contours are small, that is at most of the order log N. This is the picture of low tem- 
perature excitations of canonical ground states, and it has been treated by the authors as 
such. In particular, the entropic factor has been frequently suppressed by the microscopic 
energy cost. However, exact asymptotic results on the level of a microscopic justification 
of the Wulff construction depend, even at very low but still non-zero temperatures, on a 
non-trivial entropy/energy competition, and, hence, could not be derived in this way. 

Then there followed 15-20 years of a relative stagnation, the only contributions to the 
area being confined to generalizations of [MSI, MS2] to more complicated models [Kulj. 
A popular interest to the problem has been revived towards mid-eighties in the framework 
of an on-going mingle between probability and statistical mechanics [3c], [FO], [LeSc], 

crsT) . 

A breakthrough occurred around 1989, when Dobrushin, Kotecky and Shlosman found 
a way to derive the Wulff shape in a scaling limit of the low temperature 2D Ising model. 
They found much more: Essentially the monograph DKSfi sets up a comprehensive math- 
ematical theory of phase segregation. This theory happened to be an intrinsically prob- 
abilistic one. The DKS approach is, above all, to quantify the phenomenon of phase 
separation in terms of probabilistic limit theorems and, accordingly, to study the proba- 
bilistic structures related to the canonical states. Thus, in a sharp contrast with most of 
the preceding works, the ideology of [ DKS| ] has been from the start a very robust one and, 
actually, pertained to the whole of the phase transition region. It could be implemented, 
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however, only at very low temperatures, since the authors used low temperature cluster 
expansions as the principal tool for proving the corresponding probabilistic theorems. 

The ideas of [DKS| did not wait long to inspire a wave of investigations, even before the 
draft of the work started to circulate. Two subsequent works of a fundamental importance 
[Pf2], where an alternative simplified proof of parts of the DKS results has been 
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given using techniques, which are specific to the 2D Ising model, like self-duality, and 
[ACC|, where the Wulff construction has been derived in the context of the 2D Bernoulli 
percolation, but in a completely non-perturbative fashion, that is down to the percolation 
threshold 1/2. In both instances the exact canonical setting has been substituted by shell- 
type integral constraints, and, respectively, softer integral type limit results have been 
used instead of the local estimates of the original DKS theory. 

The results and techniques of [ACC| and [Pf2] have been combined with profound 
renormalization ideas of | Pil | and lead to an extension of this weak integral approach 
to the Wulff construction in the whole of the 2D Ising phase coexistence limit [£l]] , Q . 
Simpler proofs of some of the basic estimates of these two works (e.g estimates in the 
phases of small contours or skeleton lower bounds) have been found in [jScSl ] , [t5cS2|| , and 
the integral version of the two-dimensional DKS theory has been essentially completed 
in | P V2 1 , the estimates of the latter work being already optimal along the lines of the 
integral approach. Furthermore, Pfister and Velenik | PV1| , PV2] investigated the effect 
of boundary conditions, and in particular studied the effect of an arbitrary boundary 
magnetic field, thus providing a derivation of the Winterbottom construction. 

In spite of these successes, a non-perturbative treatment of the full DKS theory was still 
out of reach, because a key ingredient was missing: only rough estimates were available 
in the phase of small contours. By proving a local limit theorem in the phase of small 
contours, Ioffe and Schonmann were finally able to provide a non-perturbative version of 
the strong Wulff theory JESc| | . The techniques of [ISc] are based on improved versions of 
asymptotic expansions in metastable cutoff phases developed in [ ScS3| ] . 

In principle, the two-dimensional DKS theory should lead to exact expansions of canon- 
ical partition functions up to zero-order terms. This, however, requires a superb control 
over the statistical behavior of microscopic phase boundaries, which is currently beyond 
the reach for the Ising model at moderately low temperatures. A certain progress, though, 
has been reported at very low temperatures [DH], [|h| or either in the case of simplified 
models [HI]. Finally, it should be noted that at moderately low temperatures the success 



of the DKS theory in two dimensions has been by and large confined to the Ising and per- 
colation models, and that there are serious technical and possibly theoretical challenges 
to extend it to more general two-dimensional models (see Section [O] for more on this). 

On the other hand, as it has been communicated to us, an appropriate version of the low 
temperature DKS theory (as originally developed in ] DKS| ), should apply to any 2-phase 
model in the realm of the Pirogov-Sinai theory |Sh|. 

There is a strong interplay between dynamical properties of the Ising model and its 
behavior in equilibrium : in absence of phase transition, the correlations at equilibrium 
are related to the exponential relaxation of the system; instead as a phase transition 
occurs, the dynamics is driven by the evolution of droplets (nucleation, motion by mean 
curvature ...). We will not enter into details and simply refer to the seminal paper on 
metastability by Schonmann and Shlosman [ ScS3| and to the lecture notes by Martinelli 
|Ma{ (and references therein) for a survey of the recent works. Let us just mention that, as 
far as phase coexistence is considered, many dynamical results are only valid in dimension 
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2 because of the absence of a precise description of the equilibrium properties in higher 
dimensions. 



If the 2D case was subject to rapid progress, the best results for higher dimensions 
remained for a long time those of Minlos and Sinai. 

The turning point of the latest developments should be traced back to the seminal 



works by Pisztora [Pil] and by Cassandro and Presutti [CP], where crucial renormalization 
decoupling estimates have been established in the case of the nearest neighbour Ising and, 
respectively, Kac interactions. 

The basic philosophy of the Li-approach has been originally developed in the works 
[ ABCPfl , iBCPj ], RBBBPfl , |BBPf| in the context of the Ising systems with Kac pot entials , 
and, in a less explicit way, elements and ideas of the theory already appeared in [ ACCfl , 
[ Pil , || and (PV|. 

Using an embedding of the renormalized observables into a continuum setting, Alberti, 
Bellettini, Cassandro and Presutti [ ABCP| , BCPj [ emphasized the appropriateness of geo- 
metric measure theory setting, introduced relevant analytic approximation procedures (see 
Subsection 2.6.1) and proved large deviation bounds for the appearance of a droplet of 
the minority phase in a scaling limit when the size of the domain diverges not much faster 
than the range of the Kac potentials. In this scaling the system can be controlled by a con- 
tinuum limit via the T-convergence of functionals associated to the spins system JABCP ] 
and by compactness arguments [ BCP |. 

The approach of [ ABCP] and [BCP| has been extended by Benois, Bodineau, Butta and 
Presutti [BBBP], [BBP] to the case when the range of the interaction remains fixed and 
does not change with the size of the system. The latter works are, already, structured in 



a way very similar to the one we expose here. Thus the main steps of [BBBP[ and [BBP] 
comprise the coarse-graining of the rescaled magnetization profiles by the Li-proximity 
to various continuum sets of finite perimeter, surgery procedures to confine interfaces to 
tubes around the boundaries of such sets and exponential tightness arguments to reduce the 
combinatorial complexity of the rescaled problem. The essential model-related input has 



been provided by the decoupling estimates on the renormalized magnetization CP], | BZ | 
and by the result on the instanton structure of Kac interfaces DOPT1 , DOPT2]. The 
latter structure, however, yields only approximate bounds at each fixed finite interaction 
range. Consequently, the exact (van der Waals) surface tension could be recovered only 
when the range of the interaction tends to infinity, that is only in the Lebowitz-Penrose 
limit. Nevertheless, at long but finite range interactions one could say that the typical 
mesoscopic configurations concentrate on droplets with Li-almost spherical shapes. 

A complete picture of the higher-dimensional Lj-Wulff construction has been, for the 
first time, grasped and worked out in a recent remarkable work [Cel], where the corre- 
sponding results have been established in the context of the super-critical 3-dimensional 
Bernoulli bond percolation. Using novel and unusual renormalization procedures based 
on the decoupling results of [ Pil | , he has essentially rediscovered all the main steps of the 
Li-approach as described above. The main turning point of [Cel] was the introduction of 
an alternative ingenious definition of the surface tension which happened to be compatible 
with the setup of Li-renormalization procedures []. 



The work of [Cel] triggered a wave of new investigations. In [Bl] his ideas on how to 
define and treat the surface tension have been combined with an appropriate adjustment of 



It should be noted, though, that despite relative technical simplicity of this observation, the work Cel 
most certainly prompted the completion of the Li-theory by many years. 
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the renormalization approach of [BBBP] and [ BBP| ], which lead to a relatively short proof 
of the Li-Wulff construction for the nearest neighbour Ising model in three and higher 
dimensions and at sufficiently low temperatures. Most recently, a similar construction has 
been established up to the FK slab percolation threshold in [CePi|. In the latter article 



new and important techniques have been developed in order to go around mixed boundary 
conditions via bulk relaxation properties of the FK-measures. 

Although the techniques of the Li-theory might look "soft" when compared to the 
local limit setting of the DKS approach, one should bear in mind that there is always a 
"hard" step needed to initialize the Li-machinery: The renormalized mesoscopic phase 
labels have to possess sufficiently good decoupling properties. For the case of Kac models 
the corresponding estimates have been established in [CP|, [BZ], [BMP], and in the case 
of percolation (including FK for the nearest neighbor Ising model) models in dimension 
d ^ 3 in ||Pil|| , on which both [Cel],[CcPi] and [Bl] rely in a fundamental way. 

Higher dimensional Winterbottom type shapes have been recovered in the context of 
effective interface models [BI 



[DG1], 



following the original two-dimensional 



model defined and studied in |CDR|. 

The results of these works have been also formulated in terms of Li concentration 
properties, but the corresponding approach is quite different from the one we expose here. 
Thus, the analysis of | pi| | heavily relies on specific properties of Gaussian interactions. 
It should be noted, though, that, unlike in the nearest neighbour higher dimensional 
Ising case, there is better insight into the fluctuation and relaxation properties of higher 
dimensional microscopic interfaces 



_^5[, [DG1]. On the other hand, the shapes produced 
by the effective interface models are much less "physical", in particular the equilibrium 
shapes are not scale invariant, and the corresponding surface tension is not convex. 
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Part 2. Li-Theory 

On the macroscopic level the phenomenon of phase segregation is studied in terms of 
concentration properties of the locally averaged magnetization. Statistical properties of 
the microscopic phase boundaries are waved out, and the backbone of the Li-theory are 
hard model-oriented renormalization estimates, which enable a sharp surface order analysis 
of the mesoscopic magnetization profiles. Example of such coarse graining procedures in 
the case of Kac, percolation and Ising models are given in Section [2.3| . 
The averaging is performed on various mesoscopic scales: 

Mesoscopic Notation. All the intermediate scales are of the form 2 k ,k G N. For any 
M = 2 k fixed we split the unit torus T d into the disjoint union of the corresponding 
mesoscopic boxes, 

f d = \/ M k (x), (2.0.1) 



where is the scaled embedding of the discrete torus Tm = {1, . . . , M} into T as 

and, given x G T d the box M k (x) C T d is defined via 

B fc (x) = x + 



2fc+i ' 2 k+l 



Let us use J- k to denote the (finite) algebra of the subsets of T d generated by the parti 



tion (2.0.1) . Given the size of the system N = 2 n , the local magnetization Ai k on the 
M = 2 k ^ N scale is always an ^-^-measurable function. This notation should not be 
confusing: the subindex k in M k measures the "coarseness" of the mesoscopic magnetiza- 
tion profile. Thus, Mq corresponds to the microscopic configuration, and M. n identically 
equals to the averaged total magnetization. In general the local magnetization M k is a 
piecewise constant function on T d defined as 

Vxel*_ k ,VyeB n - k (x), M k (a,y) = ± £ a,. 

Notice that the microscopic counterpart of the box M n _ k (x) is the box Bm(2"x) of side 
length M centered in 2 n x. 



We formulate all the results of Section |2.1| for the nearest neighbor Ising model. Along 
with the super-critical Bernoulli percolation this is the only instance when a relatively 
complete Li-theory has been developed. In both instances, the validity of the Li-Theory 



hinges in a crucial way on the validity of Pisztora's coarse graining [Pil|, which is by 
far the most profound model related fact employed. Nevertheless, the approach itself is 
rather robust, and in subsequent Subsections we shall try to distinguish between specific 
model dependent properties and more general results. In particular, compactness prop- 
erties of local magnetization profiles are discussed in Section |2.2j without any reference 
to specific models. Instead we briefly indicate how the conditions of the corresponding 
general exponential tightness Theorem could be verified in several particular cases. 
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2.1. Results and the strategy of the proof 

2.1.1. Main results. For simplicity, we restrict to the case of the torus Tn and denote 
by hn the Gibbs measure with periodic boundary conditions. 
Define the total magnetization Mf^ as 



Mr -v = Nd Yl ai ~ 



Let us define also the set QS P as 



!B P = {[3 : Pisztora's coarse- graining hold for the Ising model at inverse temperature /?}. 



We refer to the original article [Pil] and JjCePj] for the precise relevant definitions (see 
also remark at the end of the Subsection |2.3.3[ ). It is known that 23 p contains all except 
for at most countably many points of the interval ]/3 c , oo[, where (3 C is the so called slab 
percolation threshold, which is conjectured to coincide with (3 C . 

A compact way to state the main result of the Li-theory is: 
Theorem 2.1.1. For any (3 £ 23 p and m in ]fh,m*[ 

jfc, ]v^T 1o S^(I M Tjv| < m ) = ->V/3(/C m ), 
where m = m((3) and JC m were defined in Subsection 



Remark. The above Theorem has been established for (3 3> 1 in [Bl]. The only additional 



ingredient required for an extension of the results of the latter paper to the whole of the 



temperature range j3 c was the validity of the Lemma 2.4.1. Such a statement happens to 



be highly non-trivial, and it has been proven in [CePi] along with an alternative derivation 
of the claim of Theorem |2.1.l| . □ 

Theorem |2.1.l| looks like a surface order large deviation principle. Such an appellation, 
however, would not help to explain the structure of the underlying phenomena. In fact 
Theorem [2.1.l| is essentially equivalent to a seemingly stronger statement on the macro- 
scopic geometry of the phase segregation of local magnetization profiles under the condi- 
tional measure [in ( • [Mt^I ^ rnj: 

For any function v in L 1 (T c( , [—^7*, the ^-neighborhood of v is denoted by V(v,5) 



V(« 1 i) = (o , eL 1 (f i ,[-ii ; ]) I / \v' x -v x \dxKs). 

The Li-Theorem on the phase separation says that for (3 large enough with fJ,N (■ \Mj N | < 
probability converging to 1, the function Aik is close to some translate of the Wulff shape 

More precisely, fix a number v < 1/d. 

Theorem 2.1.2. For any (3 £ 5S p and m in ]m,m*[ the following holds: 
For every 5 > 0, one can choose a scale ko = ko(f3,8), such that 



lim min fj, N I G I) V(l Km+x ,S) 
N— >oo k ^ k < vn \ m* w 



\ M -1n\ <m = 1, 
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where m and K. m were defined in Subsection 



The proofs of Theorems |2.1.1| and 2.1.2| are similar and are divided into 2 steps. The 



first step amounts to prove a compactness Theorem and the second one to derive precise 
logarithmic asymptotics. 



2.1.2. Exponential tightness. Recall [EG] that for any o positive, the set 



K a = {v £ BV(T d ,{±l}) | V({v = l})^a}, 
is compact with respect to convergence in L 1 (T rf ). 

Proposition 2.1.1. Let (3 be in 53 p . Then there exists a constant C{(3) > such that for 
all 5 positive one can find ko(5) 

Va > 0, limsup — ^-r max logfi N ( ^ £ V(K a ,d) c ) ^-C(f3)a, 

where V(K a ,5) is the S -neighborhood of K a in L 1 (T d , [— -K-, i]). 

This proposition tells us that only the configurations close to the compact set K a have 
a contribution which is of the surface order. This statement reduces the complexity of the 
problem : as K a is compact, it is enough to derive the leading terms in the logarithmic 
asymptotics for the probability of a finite number of events. 

In Section 2.2, we prove that the analog of Proposition |2.1.1| holds for a broad class of 
models. 

2.1.3. Precise logarithmic asymptotics. As the minimizers are known, it is sufficient 
to derive a lower bound for configurations concentrated close to JC m . 

Proposition 2.1.2. Let (3 be in 53 p and let m be in ]m,m*[ 

liminf — Lj min log m G V(l Km ,5)) > - Wp(K, m ) - o{5) , 

N^oo J\l a 1 ko{S) < k < un \m* J 

where the function o(-) depends only on (5 and vanishes as 5 goes to 0. 

According to proposition [2.1.1| , we will prove the upper bound only for a restricted class 
of events 

Proposition 2.1.3. Let (3 be in 53 p . Then for all v in BV(f d ,{±l}) such that Wp(v) is 
finite, one can choose 5q = 5q(v), such that uniformly in 5 < 5q 



lirnsup -n4rr . ,„ m , ax , lo S A*JV G V(v, S)j < -Wp(v)+o(5) 



1 

"nZ^ N 4 - 1 ko(5) <T<i/n 



where the function o(-) depends only on (3 and v and vanishes as 5 goes to 0. 

The Propositions above ensure that given a precision 5, there is a finite scale Aft (5) after 
which the phases are uniformly segregated with this precision. 
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2.1.4. Scheme of the proof. The scheme of the proof is well known in the soft context 
of large deviations: one first proves an exponential tightness property and then a weak 
large deviation principle (Proposition 2.1.2] holds also for any bounded variation function 
with finite perimeter). To be sure, the proof itself has nothing to do with the theory of 
large deviations: the central tools here are the renormalization estimates leading to Peierls 
type bounds and estimate in the phase of small contours, and, of course, the identification 
methods to produce the macroscopic surface tension in the precise logarithmic asymptotics. 

Thus, Proposition |2.1.1| tells us that, under the appropriate renormalization, the occur- 
rence of many small contours or of very large contours is unlikely. It is a straightforward 



consequence of the general exponential tightness Theorem 2.2.1, which we state in Sec 



tion 2.2. The statement is reminiscent to the results proven in [BBP], but the proof itself 
is based on the analysis of the phase of small contours developed in |l2| , [3cSl], [jPV2| . 

To prove Propositions [2.1.2 and [2.1.3 , we first consider the macroscopic event {^r S 
V(v, <5)j and by using several localization procedures, we reduce to compute the probability 
of microscopic events from which, adopting the procedure developed in [Gel], we can derive 
the exact surface tension factor. This enables us to avoid the computations related to the 
microscopic phase boundaries at, however, a principal cost of loosing track of the latter. 

Since the most likely configurations in {^r £ V(i>, <5)} are those for which both phases 
coexist along the boundary of d*v, we would like to prove that a microscopic interface is 
localized close to the boundary. To derive the lower bound (Proposition p. 1.2 ), one can 
enforce such a microscopic interface and then recover the surface tension factor. 



This is not the case for the upper bound (Proposition 2.1.3) because the Li constraint 
{ £ V(v, <5)} imposed on the magnetization is not strong enough to localize the interface 
close to d*v : there might be mesoscopic fingers of one phase percolating into the other. 
To circumvent this problem, we follow an argument developed in | BBBP| and first prove 
a weak localization on a mesoscopic level. This involves a surgery procedure called the 
minimal section argument. This procedure ensures that one can chop off the mesoscopic 
fingers without changing too much the probability of the event and therefore localize the 
interface on a mesoscopic level. The renormalization is an essential feature of this proof. 
Once the interface is localized on the mesoscopic level, it remains to identify surface 
tension. 



We now proceed by first defining a coarse graining and deducing the exponential tight- 



ness from Theorem 2.2.1. Then we compute the logarithmic asymptotics. 



2.2. Coarse graining and mesoscopic phase labels 

At every mesoscopic scale M = 2 k the local magnetization M.^ gives a coarse grained 
representation of the system. Statistical properties of the microscopic configurations are 
waved out, and instead one keeps track only of the local order parameters over the cor- 
responding mesoscopic blocks. These are quantified by three values ±1 and according 
to whether they are sufficiently close to one of the two equilibrium values ±m* or not. 
0-blocks play the role of the mesoscopic phase boundaries, and the ±1 blocks of the cor- 
responding mesoscopic phase regions. Thus, the outcome of the renormalization could be 
schematically represented as the following two-step diagram : 

C Microscopic f Local "| f Mesoscopic "| 

| configurations J 1 magnetization J [phase labels/ 
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There are two principal results to be discussed in this Subsection: we show that the Li- 
difference between the local magnetization and the corresponding phase labels vanishes 
on the exponential scale, and we give a general exponential tightness criterion for families 



of {±1, 0}-valued phase label functions. In Section 2.3, we will indicate how to construct 



phase labels in the case of Kac, percolation and nearest neighbor Ising models. 

Definition : A {±1, 0}-valued function u on T d is called a mesoscopic phase label, if 
there exists fcsN, such that u is an .T^-measurable function. 



2.2.1. Tightness theorem for mesoscopic phase labels. We fix now a sequence of 
non-negative numbers {p k } such that 

lim p k = 0. (2-2.1) 

k^oo 

The following compactness result holds uniformly in the microscopic scales N = 2 n . 

Theorem 2.2.1 (Tightness of Mesoscopic Phase Labels). Let N = 2 n and assume that 
{uk(oj, x)} is a sequence of random mesoscopic phase label functions defined on the common 
probability space (Qn,An,'$ > n)) such that the realizations of Uk £ J~n—ki k = l,...,n, and 
for every k the following two conditions hold: 

A. The distribution of the family of random variables {\u k {oj , x)\} j d is stochastically 

n—k 

dominated by the Bernoulli site percolation measure IPpcrc on T^_^,. In particular, 

F N (u k { Xl ) = 0,...,u k {x e ) = 0) ^(p k ) e . (2.2.2) 

B. If for two different points x,y S T^_ fc the corresponding u k -phase labels have opposite 
signs, that is if u k (x)u k (y) = —1, then on any finer scale k' ^ k any ^-connected chain of 
B n _fc/ blocks joining M n _ k (x) to M n _ k (y) contains at least one block with zero k' -label. 

Then for every a > and 5 > there exists a finite scale ko = ko(5), such that 

1 f n x<yn—dko ~] 

^logP^vKGV^^r) ^ -ci(d)min|^\ ^n^,^^} , (2-2.3) 

for all k ^ ko . 

Remark . The proof of this general theorem is given in Appendix A. Notice that for 
N sufficiently large we obtain a simpler surface order estimate which, for every v < 1/d 
fixed, holds uniformly in all mesoscopic scales ko(5) ^ k ^ vlogN, 

j^logF N (u k eV(K a ,2Sr) ^-ci(<i)^-. (2.2.4) 

Also an inspection of the proof shows that the tightness of the phase labels on a certain 
scale k does not depend on the validity of Assumptions A and B on the successive scales 
k' > k. In particular, the estimate ( |2.2.4 ) is valid on fixed (large) finite scales k = ko, once 



the Assumption A is satisfied, and once any *-connected sign changing chain of /cg-blocks 



necessarily contains a 0-block. This simplified version of Theorem 2.2.1 is used in the case 



of Kac potentials which we discuss in Subsection 2.3.1. □ 
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2.2.2. Relation to magnetization profiles. The original Gibbs measure is related to 
the above abstract setting in the following way: For every N = 2 n , one constructs a 
(possibly enlarged) probability space (fijv, An, P.iv), on which both the spin variables 
cr G {-1,+1} Tjv and various indexed families {uj^} of mesoscopic phase labels are defined. 
Such construction should enjoy the following set of properties: 

CI. The marginal distribution of a under P^v is precisely [in- 

C2. For every £ > the family {itjjj} of mesoscopic phase labels satisfies Assumption A 
of Theor em |2.2.1| with the corresponding sequence {pk,(} of site percolation probabilities 
obeying (|2.2.1| ). 

C3. For every k € {0, ...,n} and £ > the local magnetization profile and the phase 
label u* k are related as follows: Pjv-a.s., 



Mk{x) - m*ul(x) 



whenever |u£(a;) 



1. 



(2.2.5) 



Notice that both functions above are .F n _fc-measurable, that is ( 2.2.5]) should be verified 



over the mesoscopic boxes indexed by the points x £ T^_ 



ifc- 



Under conditions C1-C3, given any 5 > one can choose the accuracy £ of the coarse 
graining, a finite scale ko = ko(5,j3) and a sequence of mesoscopic phase labels {it£}, such 
that for every u < 1/d fixed, 



— L-logP.vf max ||M fc -m*t4||i>(S) < - c 2 2^ d ^ n . 



(2.2.6) 



Notice that ( 2.2.6| ) holds uniformly in the size of the system N = 2 n , once Assumptions 
C1-C3 do so. 

Let us check (|2.2.6| ). By the very construction, 



\\M k -m*ul\\i < C + 



E 1 «^,(^)=o' 



Consequently, using the domination by the Bernoulli site percolation (Assumption A), 



W N (\\M k -m*u C k U>S) 



( 



IT 



<y-C 



n— fcl 



x&T d , 

n— fc 



4(x)=o > 2 



perc 




y i 



5-C 



i(x)=o > 2 



/ 



^exp{- Cl 2^- fc )log^). 



The latter estimate is of the super-surface order once p^^ <C (5 — 0/2 and k < n/d. 



2.3. Examples of mesoscopic phase labels 



We show that mesoscopic phase labels can be constructed in the case of Kac, percolation 
and Ising models. 
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2.3.1. Kac potentials. For this model mesoscopic phase labels are denned on the original 
space of spins a G {— : the coarse graining is obtained by averaging locally the 
magnetization. Recall that we are using dyadic length scales N = 2 n . 

Phase labels are constructed in three steps. First, for any integer k and £ > 0, we 
introduce the block spin variables u* k which label the boxes B n _fc according to the averaged 
magnetization over the boxes of the linear size M = 2 k . These u k are constant on each of 
the blocks B n _fc(x) with x G T^_ fc 



u\{(T,x) 



±1 if \m3 SieB M (2«a:) °» T m *\ < C, 

otherwise. 



In the Kac case we do not use Theorem 2.2.1 in its full generality, the object of the 



coarse graining is to choose a finite scale ko, such that the family of mesoscopic phase 
labels is exponentially tight in Li. Recall that the scaling parameter is chosen such that 
e = 2~ m with to large but fixed. Eventually finite renormalization scales ko are going 
to satisfy ko = m + a®, where ao depends on (3 and (, but not on to. The sign of the 
A;o-label over a box M n _i t0 (x) depends on a more refined information on the fluctuations of 
the magnetization inside the box : we choose another scale £o; lo = to — bo, where, as in 
the case of ao, the scale bo will eventually depend only on f3 and (, and define the family 
of modified block spins {s£ } on the fco-scale as 

fi C { ^ x) A f ±1 if 4 (cr,y) = ±1, V y g T d n _ eo nm n _ ko (x) 



I otherwise. 



Finally, we define the mesoscopic phase label functions {uj, (a,x)}. If ul (a,x) = 0, we 



set u k (cr, x) = 0. If x,y G ^n~k are *-neighbors, but the corresponding modified blocks 
spins satisfy (a, x) u ko (a, y) < then n^ o (cr, x) = u ko (a, y) = 0. Otherwise, we set 

A consequence of the Peierls estimate proven in [ |CP| and [ pZj | is that assumption A is 
satisfied, namely 

Theorem 2.3.1. For any (3 > 1, there exists Co = Co(P) > 0> such that the following 
holds: For any C < Co one can choose So = £o(C)j a o = a o(C) ana " ^o = ^o(C)? such that 
uniformly in the interaction parameters e = 2~ m < eq, 

Me,iV (ul ( x i) = 0,... ,ui (x r ) = o) < exp (-prj , 



where, for every fixed e = 2 m < eq, the mesoscopic phase labels are constructed 
the scales ko = m + ao(C) anc ^ lo = m ~ ^o(C)- 



0/). 



Remark. A more refined statement implying exponential decay of correlations was proven 
in [BMP]. Notice that conditions C1-C3 of the previous Section are satisfied by definition 



of the mesoscopic phase label functions. Notice also that assumption B of Theorem 2.2.1 
is automatically satisfied on the fco-scale. Thus, the family {a^ o } is exponentially tight in 
Li. "" □ 



A similar renormalization procedure was carried out by Lebowitz, Mazel and Presutti 
|LMP| | for a system of point particles in M. d interacting with Kac potentials. In this case the 
study of phase transition in the continuum is much more involved. Beyond a proof of the 
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liquid-vapor phase transition, their results provide an accurate description of the system 
in terms of mesoscopic phase labels which represent the liquid and the gaseous phases. 
Such a coarse graining should be helpful to obtain further results on phase coexistence in 
the continuum. 



2.3.2. Bernoulli bond percolation. Bernoulli bond percolation exhibits features similar 
to the Ising model as phase transition and surface order behavior in a regime of phases 
coexistence. Nevertheless, as the setting is different from the Ising model, we briefly recall 
some notation. The set of edges is E = {{x,y} \ x ~ y}, where x ~ y means that the 
vertices are nearest neighbors. An edge b in E is open if u)b = 1 and closed otherwise. 
To any subset A <g Z d , we associate [A] e the set of edges in A. The space of bonds 
configurations in A is 

Q A = {0,1} I A k For a given p in [0,1], we define the Bernoulli bond 
percolation measure on Qa by 



be[A] e 



For simplicity <J>^ denotes the measure on Qjy — ■ 

Let w be a configuration in 17, an open path (x\, . . . ,x n ) is a finite sequence of distinct 
nearest neighbors x\, ... ,x n such that on each edge coi XijX \ = 1. We write {^4 «-> B} 
for the event such that there exists an open path joining a site of A to one of B. The 
connected components of the set of open edges of u are called w-clusters. 

A phase transition is characterized by the occurrence of an infinite cluster. Define Q p 

by 

Q p = lim $^({0 ^ OTjv}) , (2.3.1) 

N— >oo 

then there is a critical value p c in ]0, 1[ such that for any p below p c there is no percolation 
and Q p = 0, instead for any p above p c the occurrence of an infinite cluster starting from 
has positive probability Q p . In the thermodynamic limit, there exists only one limiting 
Gibbs measure and almost surely a unique infinite cluster with local density Q p . In order 
to mimic the coexistence of 2 phases in the finite domains HV, we say that one phase is 
formed by the largest cluster and the other phase by the other clusters. 

For this model, Pisztora introduced a renormalization procedure [ [Pil| , |DePi| , [Pi2] 



which holds as soon as p > p c and d ^ 3. The mesoscopic phase labels {uf} will be 
defined for any mesoscopic scale M = 2 k , where k is an integer which eventually depends 
on N. This construction requires 2 steps. The first step is to retain only the main features 
of the typical configurations on finite size boxes Mm- Then we attribute a sign to the 
blocks B ra _fc according to the phase they represent. Set M' = 2M. For any x in T^_ fc , the 
following events depend only on configurations in the box HH>M'(2 n x). 

U x = G Sljv | there is a unique crossing cluster C* in M M '(2 n x)}. 

A crossing cluster is a cluster which intersects all the faces of the box. Let i be an integer 
smaller than k which will be fixed later 

Rx = U x |w G iljv | every open path in MM'(2 n x) with diameter larger than 2 e 
is contained in C* j, 
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where the diameter of a subset A of Z d is sup^g^ ||x — y\\\. Finally, we consider an event 
which imposes that the density of the crossing cluster in Bm(2"x) is close to Q p with 
accuracy £ > 

V x = U x f){u€n N \ \C* n M M (2 n x)\ G [e p - c, Q P + C] 



Each box B n _fc(x) is labeled by the variable u k (oj,x) 



where I • I denotes the number of vertices in a set. 



otherwise. 



Let {x\, . . . , x r } be vertices in T d l _ k not *-neighbors of x, then [Pil| implies that for every 
p > p c , there exists ko{p, (), and £o(p) such that for all k ^ ko and k ^ £ ^ £q 



(u<( x ) = 0\ ul( Xl ),...,ui(x r ^j < exp(- Cl 20+exp(-c 2 (C)2 fc ), 



From |LSS(| (Theorem 1.3), we deduce that for k and £ large enough, the random variables 
{u k (x)} are dominated by a Bernoulli site percolation measure wPk 



pcrc 



p k < exp(-c(C) 2*). (2.3.2) 

A straightforward way to recover the previous statement is to partition T^_ fc into e(d) 
sub-lattices (T^_ fc _ 1 J i<c ^ which are translates of T^_ fc-1 . Any collection of vertices 

{xi,...,x r } in T^_ fc can be rearrange into c(d) subsets {x^ , ...,x$} such that each 
{x\ , . . . , Xr}} belongs to T^ l _ k _ l i . Applying Holder inequality, we get 

C(d) ! 

^ N [ni(x 1 ) = 0,...,ui(x r ) = 0) < ]J& N (ui(x?) = 0,...,ui(xW) = 0) 7m . 

i=i 

As the vertices in T d l _ k _ li are not *-neighbors in T^_ fc , the domination by a Bernoulli 
product measure follows. 

We say that a block M n _k(x) is regular if u k (x) = 1. Finally we define the mesoscopic 
phase labels u k to be equal to 1 on the regular blocks connected to the largest cluster 
and to —1 on the regular blocks disjoint from the largest cluster. Otherwise, we set 
v£ k {<jj,x) = u k (u>,x) = 0. From ( 2.3.2j ), the mesoscopic phase labels satisfy assumption 
A. Notice that if x and y are *-neig hbors in T d n _ k the boxes B M ,(2 n x) and M M >{2 n y) 
overlap. Choosing the parameter £ ^ k — 3 we insure that if the boxes B n _^(x) and 
B n _fc(y) are both regular, then the crossing clusters in these boxes are connected. This 
implies that assumption B is satisfied : two blocks with A:-labels of different signs cannot 
be *-connected. 

The Bernoulli bond percolation model is precisely described by Pisztora's coarse grain- 
ing, namely on a sufficiently large scale 2 k , the typical configurations have a unique crossing 



cluster surrounded by small islands of size smaller than 2 . According to Theorem 2.2.1 
the family {u k } is exponentially tight in L 1 . 
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2.3.3. Ising nearest neighbor. An extension of the preceding renormalization procedure 
applicable to the Ising model has been also introduced in [pPil| ] . Unlike Ising model with 
Kac potentials, this coarse graining is defined on an enlarged phase space via the FK 



representation. For a review of FK measures, we refer the reader to [Pil], [A.CCN] and 

[cp 



Let us recall the definition of the random cluster measures (or FK measures) which 
are a generalization of the Bernoulli bond percolation measures with correlated bond 
distribution. To any subset A of Z d and ir included in <9A, we associate a set of edges 

[Ale = {{x,y} | a;~ V, x G A, y G AU^}, 

and the space of configurations in A is £l\ = {0, l}t A ^. The first step is to introduce a 
measure on $7^. A vertex x of A is called 7r-wired if it is connected by an open path to 
7r. We call 7r-clusters the clusters defined with respect to the boundary condition ir : a 
7r-cluster is a connected set of open edges in f2T and we identify to be the same cluster all 
the clusters which are 7r-wired, i.e. connected to tt. For a given p in [0, 1], we define the 
FK measure on Q\ with boundary conditions ir by 



*i'»=i n (i-p) i-< v*)2 

A W[A]« 



where Z^' p is a normalization factor and (^(u) is the number of clusters which are not ir- 
wired. If ir = dA then the boundary conditions are said to be wired and the corresponding 
FK measure on is denoted by < 3?^' p . Finally, the periodic measure on the torus T^r is 
denoted by <3?^ r ' p and the phase space by ft^ r . 

In order to recover the Gibbs measure we fix the percolation parameter pp = 
1 — exp(— 2(3) and generate the edges configuration u> in £}^ r according to the measure 
<F™. Given u, we equip randomly each w-cluster with a color ±1 with probability tj in- 
dependently from the others. This amounts to introducing the measure on { — 1, 1} Tjv 
such that the spin Oi has the color of the cluster attached to i. The Gibbs measure /j,n can 
be viewed as the first marginal of the coupled measure Pat(<t, u>) = P%(a)$ p *' P0 (w) on the 
space {—1, 1} Tjv (8> ti^ 1 • In the case of 7r-wired boundary conditions, the spins attached 
to the 7r-wired cluster are equal to 1. 



As a consequence of this representation, one has for any increasing sequence of sets Ajv 

m * = J im Oo) = iim $If ({o « aA^}) = e p „. 

In the following, we use m* or Q P/3 depending on the context. Furthermore, we suppose 
that 

lim $^({0 ~ dA N }) = lim c^f({0 <- dA N }) = O P0 . (2.3.3) 

This property is satisfied for all outside a subset of R which is at most countable (see 
Lebowitz and Pfister [ pflj ). 

On the scale M = 2 k , we define, in the same way as for Bernoulli bond percolation, the 
variables u^{uj,x) which are piecewise constant on each box B ra _fc(x) with x in T^_ fc . The 
mesoscopic phase labels depend on the averaged magnetization in regular blocks. Define 
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the label of 



(x) by 



u C (a,LU,x) = I s[ S n ( C *) if u C k {u,x) = 1 and \M k (a,x) 
k \ otherwise, 

where C* is the crossing cluster in B^(2 n a;). 



sign(C*)m*| < 2(, 



In a regular box B ri _fc(x) (i.e. u^{x) = 1), the averaged magnetization is controlled 
by the random coloring of the small clusters included in Mm(2 u x). So that the averaged 
magnetization in a regular box is independent of the configurations in the neighboring 
boxes. In the case of Ising model, the additional parameter I = £(k) is tuned in order to 
control the fluctuations of the magnetization over the small clusters. As a consequence 
of this, assumptions A, B and C1-C3 are satisfied for pp above a certain non-trivial slab 



percolation threshold p^ , which is conjectured to coincide with pp c (see ] Pil ] for details), 
and Theorem 2.2.1 holds. 

Remark . Using the notations of this Subsection, the set <B P introduced in Subsection 
|2.1.1| could be defined as 



«B p = {13 : (3 > p c and fl2lOD holds} 



2.4. Surface tension 



We are going to derive Propositions |2.1.2 and 2.1.3 for Ising model with nearest neighbor 
interaction. As explained before, the philosophy of the proof is to start from the macro- 
scopic level and to localize successively on finer scales with the help of a coarse graining. 
The approach itself is quite general. Nevertheless the coarse graining is model dependent, 
therefore we will need first to state an alternative representation of the surface tension in 
terms of the FK representation in order to use the estimates which will be obtained from 
Pisztora's coarse graining. The idea of such definitions has been introduced in [Cel|. 



2.4.1. FK representation. We fix n a vector in § and study Tp(n). Following no 



tation of Subsection 1.2.4 , we consider, for any e positive, the parallelepiped A(N,eN) 
of M. d oriented according to n. Namely, the basis of A(N,eN) with side lengths equal to 
N is orthogonal to ft and the other sides have lengths equal to eN. For simplicity its 
microscopic counterpart A(N,eN) n Z rf will be denoted by Ajv(e)- 

By using the correspondence between the Ising model and the FK representation, one 
can rewrite t$ in terms of the bond model. Let {d + A^(£) <A <9~Ajv(e)} be the event such 
that there is no open path inside An(e) joining d + Ajy(£) to d~AN(e). Then, 

r p {n) = lirn^ -^log^^^a+A^Ce) * (TA^e)}). (2.4.1) 



Notice that the event {<9 + Atv(£) jA d~A^(e)} takes only into account the paths inside 
A7\r(e) and not the identification produced by wired boundary conditions. The relation 
above will be useful only in the proof of Proposition |2.1,2j . 



We are now going to state an approximate expression of the surface tension which 
is weakly dependent on the boundary conditions. It will be used in the derivation of 
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Proposition 2,1.3 , Let A' Ar (e) be the the parallelepiped 

A' N (e) = {icA N (e) \ i- n € [~N, |iV]} , (2.4.2) 

and denote by d top A' N (e) (resp d hot A' N (e)) the face of d + A' N (e) (resp 9~A' Ar (e)) orthogonal 
to ft. Let {3 top A^ v (e) <A O bot A' Ar (e)} be the event such that there is no open path inside 
A' N (e) connecting d top A' N (e) to d hot A' N (e). One has 



Lemma 2.4.1. []B| /? > 1, ]CeP| (3 G *8 P ] For any € <B 



Mn) = ({d top ^ N (e) <h d hot A' N (s)}) +c^(vr), (2.4.3) 

where the function c E) n goes to as N tends to infinity and e goes to 0, uniformly over 
the boundary conditions it and n 6 S rf_1 . 

As it will be explained in Part ||] on the wetting phenomenon, the system is in fact ex- 
tremely sensitive to boundary conditions. Nevertheless in the above Lemma, the interface 
is constrained to be in A' N (e), so that it does not feel the influence of the boundary : the 
boundary conditions are screened because the system relaxes to equilibrium in the region 
A N (e)\A> N (e). 

Let us first examine the influence of the boundary conditions ir on the faces of An(e) 
orthogonal to ft. As {d to PA' N (e) <f* d bot A' N (e)} is a decreasing event, FKG inequality 
imply that it is enough to check that 

Mn) = -^ =I log$^g({5+A' JV ( £ ) t4 d~A' N (e)}), (2.4.4) 

where $j^J?J\ is the FK measure with free boundary conditions on the faces orthogonal to 
ft and wired on the others. This can be proved by means of a Peierls argument for (3 large 
enough [Bl] or by an analysis of the relaxation of the clusters density for (3 in Q3 P [CcPi]. 



As already noticed in [pel | in the context of percolation, the influence of the boundary 



conditions on the sides of An(s) parallel to ft is negligible as e goes to 0. This explains 
that the factor c £j at(-) vanishes uniformly over the boundary conditions. 



2.4.2. Extended representation. We would like to stress that the previous treatment of 
the surface tension is not satisfactory and a more coherent approach would be to consider 
a more general definition independent of the model in terms only of mesoscopic phase 
labels. In fact, a definition of surface tension valid in an abstract setting would be difficult 
to use because the surgical procedure of the minimal section argument requires a precise 
knowledge of how the microscopic system is related to the mesoscopic phase labels. 



2.5. Lower bound : Proposition 2.1.2 



The proof is divided into 3 steps. We first start by approximating the surface d*JC m by 
a regular surface dK and imposing the condition that a mesoscopic interface exists close 
to dK. Then, using the definition of surface tension ( 2.4.1] ), we derive Proposition [2.1.2 , 
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Figure 6. Polyhedral approximation. 

2.5.1. Step 1 : Approximation procedure. A polyhedral set has a boundary included 
in the union of a finite number of hyper-planes. The surface d*fC m can be approximated 
as follows (see Fig. ||) 

Theorem 2.5.1. For any 5 positive, there exists a polyhedral set K such that 
% e V(l Km ,5) and \Wp(K) - W p {K rn )\ < 8. 

For any h small enough there are I disjoint parallelepipeds R 1 , . . . ,R e with basis B 1 , . . . , B £ 
included in dK of side length h and height Sh. Furthermore, the sets B 1 , . . . , B e cover dK 
up to a set of measure less than 5 denoted by U s = dK \ Ui=i B l and they satisfy 



e 

i=i jBi 



^6, 



where the normal to B l is denoted by Hi. 

The proof is a direct application of Reshtnyak's Theorem and can be found in the paper 



of Alberti, Bellett ini []AlBo | 



Using Theorem |2.5.1| , we can reduce the proof of Proposition 2.1.2 to the computation of 
the probability of {— r € V(J^-,<5)}. According to ( |2.2.6 ) the estimates can be restated in 
terms of the mesoscop ic ph ase labels. For any 5 > 0, there exists ( = ((5) and ko = ko(S) 



such that Proposition |2.1.2j will be implied by 
1 



liminf — -j—r mm 

N->oo N"-- 1 k (8) ^ra 



logPjv(i4<E V(l£,5)) ^ 



W p {K)-o{5). (2.5.1) 



2.5.2. Step 2 : Localization of the interface. The images of K, R i and U 5 in T N will 

be denoted by Kn, R % N and Ufr. In order to enforce a mesoscopic interface which crosses 
each R l N , we define the event 



A = f]{d + R l N ^8-R l N } 



We consider also B the set of configurations such that the bonds at distance less than 10 
of are closed. Notice that these events depend only on bonds variables. One has 



N 



u£€V(l#,tf)) ^Piv({4 G V(l R ,5)}nAnB 



(2.5.2) 
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The interface imposed by the event APiB decouples Kn from its complement, therefore 
the system is in equilibrium in Kn and : a proof similar to the one of Theorem 2.2.1 
implies that one can choose £' = C'(<5) an d k' = k' (5) such that 

lim max Pjv ( f \ul (x) — 1| dx ^ - or f \ul (x) + 1| dx > - AC\ B J = , 

N-^oo k' {S) <to V^A 2 JA 2 / 

where A stands for .ff or AT C . So that ( |2.5.2| ) can be rewritten for TV" large enough as 

min F N (u{ € V(%, 5)) > ]- <I>™ (AnB). (2.5.3) 

ko{6) ^ k ^ vn V / o 



2.5.3. Step 3 : Surface tension. Combining the definition of surface tension ( |2.4,1 ) 
inequality (|2.5.3| ) and Theorem 2.5.1 , we get 



l t [ni n^i , tK ^< logFN H G V(I ^' *)) > - E /- Tp(fii) dnt 1 - o(5). 

N—kx iV" x ko(S)^k^i/n \ 1 / ~i J B i 



t=l 

d-l 



We have also used the fact that the event £> is supported by at most c(d, 5)N edges 
where c(cf, (5) vanishes as 5 goes to 0. Therefore the probability of ,6 is negligible with 
respect to a surface order. 



2.6. Upper bound : Proposition 2.1.3 



The proof is divided into 3 steps. First we decompose d*v in order to reduce the proof 
to local computations in small regions. Then in each region we localize the interface on 
the mesoscopic level via the minimal section argument. Finally the last step is devoted to 
the computation of the surface tension factor. 

2.6.1. Step 1 : Approximation procedure. We approximate d*v with a finite number 
of parallelepipeds (see Fig. |7|). 

Theorem 2.6.1. For any 5 positive, there exists h positive such that there are £ disjoint 
parallelepipeds R} , . . . ,R included in T d with basis B 1 , . . . , R of size h and height Sh. 
The basis B % divides R % in 2 parallelepipeds R l ' + and R l ~ and we denote by Hi the normal 
to B % . Furthermore, the parallelepipeds satisfy the following properties 



[ |%(ar) - v(x)\ dx < 5vd(&) and I V f Tp{iU) dH {d ~ 1] - Wp(v) 

JR? 1 ~l J Si 

where = — an d the volume of R % is vol(i?*) = 5h d . 

This Theorem is a rather standard assertion of the geometric measure Theory. A variation 
of it has been formulated and applied in the context of the Li-theory of phase segregation 



in [ABCP] along with a sketch of the proof, which, however, contained a gap (see | B1 | 



for a detailed proof along the lines of ABCP|). A very clean alternative derivation of a 



similar result has been given by Cerf [ Cel l using the Vitali covering Theorem. 



Theorem 2.6.1 enables us to decompose the boundary into regular sets (see Fig. |?]) so 



that it will be enough to consider events of the type 

h e nv ( ff,M«')) 

l i=l 
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Figure 7. Approximation by parallelepipeds, 
where V(R l ,s) is the e-neighborhood of 

V{R\e) = W E L 1 (T d ) | J \v'(x) - X^(x)\ dx ^ e 

Using (|2~2l| ), we see that to derive Proposition [2.1.3| , it is equivalent to prove the following 
statement for any 5 positive and ko = ko(5), C = C($) 

1 £ 
limsup — -r-r max log Pat (u£ G P| V(-R l , <5vol(E 1 ))) < - + C(fi,v)6. 

2.6.2. Step 2 : Minimal section argument. The microscopic domain associated to R 1 
is R% = NR { DT N . We also set R^ = NR*>+ D T N and R%~ = R%\ R^ . At the scale 
M = 2 fc , we associate to any configuration (a, u>) the set of bad boxes which are the boxes 



intersecting R b N labeled by and the ones intersecting R^ (resp R^ ) labeled by — 1 
(resp 1). For any integer j, we set B l, i = B % + j c(d)2 n ~ k Hi and define 

B% j = {f e R% | 3x € B^, ||/ - iVs[|i ^ 10}. 

Let B\ be the smallest connected set of boxes Mm intersecting B 1 ^ . By construction the 
B\ are disjoint surfaces of boxes. For j positive, let nf(j) be the number of bad boxes in 
B\ and define 

n+ = min{n+(j): <j<^2- fe }. 

Call j + the smallest location where the minimum is achieved and define the minimal 
section in it^ + as B\ . For j negative, we denote by B\ the minimal section in R l ^~ and 
the number of bad boxes in B\ (see Fig. |8|). 
For any configuration (a, u>) such that u" k {a, uj) belongs to Di=i V(R i ,S\ol(R 1 )), one can 
bound the number of bad boxes in the minimal sections by 

i 

nf + n~ ^ 5C X (v^- 1 ^^ . (2.6.1) 

i=l 

Such an estimate implies that a mesoscopic interface is mainly located between the 2 
minimal sections and that only some mesoscopic fingers attached to the interface may 
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bad blocks 




{« = !} 



bad blocks 
Figure 8. Minimal sections. 



percolate. As these fingers will cross the minimal sections through bad boxes, the strategy 
is therefore to modify the configuration tu on the bad boxes so that no fingers can percolate 
in the new configuration. More precisely, we introduce the set 



3a such that ut(a, uj\ 



and for any u> in A define u> the configuration with closed edges on the boundary of the 
bad blocks in the minimal sections and equal to u otherwise. Inequality ( |2.6.1 ) implies 
that u! and Hi differ only on at most SC2(v)N d ~ 1 edges, so that we can control precisely 
the cost of the surgical procedure which consists in isolating the bad blocks in the minimal 
sections by closing the edges around them. 



€ P| V{R l ,5vol(R 1 )) 



(A) 

exp(5C 3 (v,P)N a 



(2.6.2) 



- 1 ) *r p ^), 



where A = {u> \ u> £ A}. 



2.6.3. Step 3 : Surface tension estimates. Let R 1 ' be the parallelepiped included in 
R % with basis B % and height Its microscopic counterpart is R % N ' . We are going to check 
now that A is included in P|^ =1 {9 top R l N ' <A d hot R % N '}. This amounts to say that not only 
the minimal section argument enables us to find a mesoscopic interface in R l N but that 
in fact this interface exists on the microscopic level. To see this, choose any configuration 
u) in A which contains an open path C joining d to PR%' to d bot i4 and 

suppose that C 

crosses the minimal sections without intersecting a bad box. Then C intersects 2 regular 
boxes Bjv/(2 n x + ) and Mm (2 n x~) in B\ and B\ . According to the definition of the coarse 
graining, this would imply that the crossing clusters of Bm(2"x + ) and Bjv/(2 n x _ ) are 



connected to C, so that ifi k {x + ) 
box. 



u^ix ). Therefore one of these boxes has to be a bad 
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From we get 

P ^v f 4 e H V ^'' 5vol (^))l < exp (5 C 3 (v, l3)N d - 1 ) 



8=1 



1=1 

Conditioning outside each domain R l N and using ( |2.4.3 ), we derive 

limsup — max log Pjv (u£ G P| V(fl\ tfvolfH*))^ < 



1=1 



- r 

I Tp(n i )dH x + C 4 (p,v)5. 

8=1 



This concludes the Proposition. 

2.7. Open problems 

We would like mention some open questions related to the Li-theory 

1. Extention of the Li-theory to general finite range models and to the context of 
Pirogov-Sinai Theory. 

2. Proof of the Wulff construction for continuum models in an Li-setting. 

3. Upgrade of the concentration properties to the Hausdorff distance, based on more 
delicate versions of the minimal section argument; some results of this sort should 
appear in |BIV]. 

4. A more challenging problem would be to provide an accurate description of phase 
segregation a la DKS. In particular one should understand how to control phase 
boundaries and prove local limit results with boundary conditions which are only 
statistically pure. 
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Part 3. Dobrushin-Kotecky-Shlosman (DKS) theory in 2D 

In this part we review and explain the results on phase separation in the two-dimensional 
nearest neighbor Ising model as enforced by the canonical constraint on the magnetization 



[ DKS ], [[Sc]. The theory is built upon sharp local estimates over finite volume vessels A^v 
and on the probabilistic analysis of the random microscopic phase separation line. We 
focus here on the "free" spatial geometry of the phase segregation, that is disregarding 
the boundary effects. These effects could enter the picture in two different ways: in terms 
of the boundary conditions on dAjy and in terms of the geometry of A^y. In the former 
case the minority phase could be absorbed by part of the boundary dAjy. This and 
related phenomena are discussed in Part 4. In the second case the finite vessel Ajy might 
not be able to accommodate the corresponding optimal crystal shape. Such a geometric 
constraint is, from the point of view of the microscopic theory, merely a technical nuisance, 
though, on the macroscopic level, it might lead to formidable variational problems. We 
go around this domain geometry issue by choosing A at to be of the Wulff shape itself 

A N = AT/dnZ 2 , 

where K,\ is the unit area Wulff shape. Thus, A^v accommodates any optimal shape of 
area smaller than N 2 . 

The corresponding finite volume canonical Gibbs measure is then defined by 



l Nt _ ( • | M N (a) = -N 2 m* + a N ) , (3.0.3) 



where Mn = J2ieA a i ^ s the total spin, m* = m*((3) is the spontaneous magnetization, 
and oat points inside the phase transition region, a at £ (0, 2N 2 m*). In the sequel we shall 
use the shortcut ^ N _ for the finite volume measure /J^ _ . 

Notation. The values of positive constants ci, C2, ... are updated with each subsection. 



3.1. Main Result 

DKS theory gives a comprehensive solution to the following problem of phase separation: 

Problem 1. For > (3 C and aj\r £ (0, 2N 2 m*) characterize typical spin configurations a 
under the canonical measure ( [3.0.3D . 

An ostensibly simpler problem is 

Problem 2. For > (3 C and £ (0, 2N 2 m*) find sharp local asymptotics of 

l4f,-i M N = -m*N 2 + a N ). 

In fact both problems are equivalent. In particular, the phenomenon behind the shift of 
the magnetization is inside the phase transition region not a bulk one (and hence is not 
in the realm of the usual theory of large deviations), and the crucial role is played by the 
spatial geometry of symmetry breaking. 

3.1.1. Heuristics. Under the finite volume pure state [iP N _ the typical maximal size of 
± contours is of order log N . One could then visualize a typical microscopic configuration 
a on Aat in terms of an archipelago of small (that is of the maximal size ~ log A) "+" 
islands which could contain still smaller "— " lakes etc. This archipelago spreads out 
uniformly over A^r, and the density of the plus "soil", which spells out in terms of the 
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magnetization M^(a) as (|Ajv| + Mjv((t)) /2|Ajy|, is close to its equilibrium value 



|Ajv| + {M N f N _ i 



m 



2\A N \ 2 

Thus, one could think of two different competing patterns behind the aAr-shifts, 0, 
of the magnetization Mjy from its equilibrium value (Mjqy N _ ~ —m*\A^\: 

1) The density of the archipelago increases in a spatially homogeneous fashion without, 
however, altering the typical sizes of the islands. 

2) Spatial symmetry is broken, and an abnormally huge island of the "+" phase of excess 
area ~ /2m* appears. 

Heuristically, the first scenario corresponds to Gaussian fluctuations, and its price, in 
terms of probability, should be of order 

exp(- Cl (f3)a 2 N /N 2 ). 

Phase segregation manifests itself in the second scenario, and the probabilistic price for 
creating such a huge island is proportional to the length of its boundary 

exp (-c 2 (/3)- v /ojv) • 

A comparison between the two expressions above suggests that the first scenario should 
be preferred whenever a at <C A^ 4 / 3 , whereas large shifts a at 3> A r4 / 3 should result in the 
phase segregation picture described in the second scenario. This indeed happens to be the 
case, and we refer to [DS| and | ISc ] for a complete rigorous treatment^. 



For the sake of the exposition, we shall stick here to the possibly most interesting case 
of a at ~ A^ 2 , which corresponds also to the macroscopic type of scaling discussed in Part 2. 
The DKS theory gives then the following sharp characterization of the phase segregation in 
the canonical ensemble: under fjP N _ ( • \Mn = —m*N 2 + a at) a typical spin configuration 
a contains exactly one abnormally large contour 7 which decouples between the "+" phase 
(inside 7) and the "— " phase (outside 7). In particular, the average magnetization inside 
(respectively outside) 7 is close to m* (respectively —m*), and the area encircled by 7 can 
be thus recovered from the canonical constraint, 

m ■ ' 2 ~ ~*- 2 ' ' ^' ~ aN 



|int (7)1 — m* (N — |int (7)]) « -m*N 2 + a N =^ |int (-7) | 



2m* 

Under the scaling of A at by l/N, that is into the normalized continuous shape /C C M 2 , 
the microscopic phase boundary 7 sharply concentrates around a shift of the Wulff shape 
of the corresponding scaled area /2m* N 2 (Fig. |9|). 

3.1.2. DKS theorem. More precisely, for any r £ M+ let K r to denote the Wulff shape 
of the area r. Also given a number s G R + , let us say that a microscopic contour 7 is 
s-large, if dianioo^) > s. 



Theorem 3.1.1 ( jDKS | ^ pSdJ ). Let the inverse temperature (3 > (3 C be fixed, and let the 



sequence {a^}, —m * N~ + £ Range(MAr), be such that the limit 

a = lim ^ G (0,2m*(P)) 



7 The critical case of ajv ~ N 4/l3 is still an open problem. 

8 In the original monograph [DKS] the corresponding results has been derived in the context of the Ising 
model with periodic boundary condition. 
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Figure 9. DKS picture under the 1/N scaling: On the left the microscopic Ajv box with 
the unique K log iV-large contour 7. On the right the continuous box /Ci with the scaled 
image of 7. 

exists. Then, 

log f4f_(M N = -m*N 2 +a N ) = -Wp(dK^j(l + O^ 1 / 2 logiV)). 

Moreover, if K = K{f3) is large enough, with fjP N _ ( • \Mn = —N 2 m* + on) -probability 
converging to 1 as N — ► 00: 

1. There is exactly one K '(/?) log N -large contour^. 

2. This 7 satisfies 

min-Lfef 1 ,x + dK ± M_) < afflN-^y/logN (3.1.1) 

X j\ 2m* 

and 

rnin-^Area^ int (7) A (x + u) ) < c 2 (/3)iV~ 3/ Vlog ^- (3-1-2) 

3.1.3. DKS theory. The DKS theory views the production of the event {Mn — m*N 2 + 
ajy} in terms of a two-step procedure: On the first stage a length scale s = s(N) is chosen, 
and all the microscopic s-large contours (7i,...,7n) are fixed. If the total area inside 
these s(iV)-large contours is smaller than /2m* , then the total magnetization still 
has to be steered towards the imposed value Mn = —m*N 2 + a^, but already under 
the constraint that all the ± contours different from (71,..., 7 n ) are s(iV)-smaIl. The 
probability A%- {^n = —m*N 2 + ajv) reflects the price of the optimal strategy along 
these lines. 

We record the two steps of the DKS theory as follows: 

1) Study the statistics of s(A r )-large contours under IjP n _■ 

2) Give local limit estimates on the magnetization in the s(iV)-restricted phases. 

The introduction of s(iV)-cutoffs leads to the separation of the length scales which has 
a double impact on the problem: it sets up the stage for the renormalization analysis of 
microscopic phase boundaries, and it improves the control over the bulk magnetization 
inside the corresponding microscopic phase regions. Let us try to explain this in more 
details: As far as the statistics of the s(A r )-large contours is considered, we are interested 
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Figure 10. Two microscopic contours 71 and 72 are compatible with the same skeleton 
S = (ui, ...u„). 



in giving sharp estimates on the [/L _-probability of the events of the type 

{ s(iV)-large contours of a encircle a certain prescribed area } . 

The point is that the contribution of any particular microscopic contour to the probability 
of such an event is negligible. In other words, one also has to take into account the entropy 
(number) of all the contributing contours. The required entropy cancelation (and hence the 
production of the relevant limiting thermodynamic quantity - surface tension) is achieved 
by means of a certain coarse graining procedure, the so called skeleton calculus, which 



we describe in Section 3.4. Roughly, instead of studying the probabilities of individual 
microscopic contours one considers the packets of all contours passing through the vertices 
of a given "s(A r )-skeleton" S = (ui,U2, ■■-,u n ) and staying within a distance of the order 
s(N) from the closed polygonal line Pol(S') (Fig. [To| ) . The distance between successive 
vertices of S complies with the length scale s(N), — «i||oo ~ S (N)- Surface tension 

is produced on the level of skeletons. In fact, the probability of observing a ± contour 
compatible with a given skeleton S admits an asymptotic with s(N) f 00 description 

4_(S) x exp{-W^(Pol(5))}. (3.1.3) 



We quote the precise result in Section |3.4| , which we devote to a general exposition of the 
skeleton calculus. 

Since the vertices of S are s(iV)-apart, and the surface tension tr is strictly positive for 
all P > p c , the energy Wp (Pol(5)) controls the number #(£) of vertices of S as 

m K C3(ffl^«. (3.1.4) 

When combined with fl3,1.3|) this leads to the reduction of the combinatorial complexity 
of the problem: the number of different skeletons of a fixed energy Wjv does not compete 
with the approximate probability exp{— WV} to observe any such skeletons. Thus, the 
study of {Mj\r = —m*N 2 + a^} reduces, in terms of skeletons, to the maximal term 
estimation. It should be stressed, however, that unlike the coarse graining procedures of 
the Li theory, the mesoscopic objects (skeletons) of the DKS theory closely follow the 
microscopic structure of phase boundaries. 
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The local limit estimates in the s(A)-restricted phases are, therefore, required uniformly 
over finite lattice domains whose boundaries are carved with s(A)-large contours compat- 
ible with not too costly skeletons. This imposes a natural restriction on the length of 



these boundaries, and we shall describe the appropriate family of domains in Section 3.2 
along with the exposition of the corresponding uniform local limit results. Intuitively, long 
contours are responsible for long range dependencies between spins, and, therefore, the 
s(A)-cutoff constraint improves the mixing properties of the system and helps to extend 



the validity of classical (Gaussian) behavior of moderate deviations. In Section 3.3 we 



quote the corresponding relaxation and decay properties which lie in the heart of the local 



limit estimates. In Section 3.5 we give an outline of the proof of the DKS theorem. 



Finally, the (long) list of open problems is briefly addressed in Section 3.6. 

3.2. Estimates in the phases of small contours 

As it has been mentioned, the estimates in the phase of small contours should be derived 
uniformly over a family of lattice domains whose boundaries are composed of not too costly 
s(A)-large contours. 

Definition Basic family T>n of subsets A C An' We fix two numbers a (small) and R 
(big). 

AeV N <=^ aN 2 < \A\ and \dA\ < RN log A. 

□ 

We fix a basic scale s(N) = A' log A of large contours, where A = K(J3) is a sufficiently 
large number, so that A log A-contours are highly improbable under the pure state fj, N _. 
Of course, exactly the same number A appears in the statement of Theorem |3.1.1| . The 
upper bound on dA in the definition of the family T>n states that the configurations with 
total length of A log A large contour exceeding A A log A are ruled out. This conclusion 



is explained in more detail in Section |3.4| (see the remark following Lemma 3.4.1 ). 



3.2.1. Structure of local limit estimates. Let us turn now to the structure of local 
limit estimates in the s(A)-restricted phases. First of all, given any A C Z 2 , the s- 
restricted phase on A is defined via 



3,8 t \ A 3 
_ ( ' ) = It a. 



All ± contours are s-small 



We would like to study the probabilities of deviations ^ of the total magnetization 
Ma from the corresponding averaged value (M^)^' s _. Let us define the set of feasible 
values of such deviations as 

M+ = {ajv^O: (M A fx"- + a N G Range(M A )} . 

Roughly, the cutoff s extends the validity of Gaussian moderate deviations for the following 
reason: The price of shifting the magnetization by a^r on the expense of s(A)-small 
contours is of the order (cln/s 2 )s ~ cin/s(N). This should be tested against the Gaussian 
moderate deviation exponent of the order a 2 N /N 2 . Thus the Gaussian behavior should 
prevail once a at <C A 2 /s(A). Of course, the latter constraint on becomes less stringent 
as s(A) decreases. On the rigorous mathematical part the classical approach to estimating 

4' s _ (m a = (M A fx s „ + aN 



amounts to first finding the value of magnetic field 

g = g(A,s(N),a N ) 
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such that the expected magnetization under the g-tilted state is precisely what we want, 

{Ma) P A%, 9 = (M A f/- + "at, (3.2.1) 
and, then, to rewrite the ^'^-probability in terms of the \i A _ g one: 

a ( - * , , , , fl.* . \ 
Ma,- 




log(e^)5 s _} AU 9 {M A = (M A \^ S 



h dhdr } ( M A = {M A )X_ >g 



(3.2.2) 



One then tries to derive sufficiently precise estimates on the semi-invariants of fJ> A _ h and 
to prove a local CLT under n A ' _ . Thus, it is extremely important to understand how the 
magnetization (M A ) A ,S _ g and other semi-invariants of ^ A s _ g change with the magnetic 
field g in the phase of s(iV)-small contours. 

Breaking of the classical limit behavior in the s(iV)-restricted phase manifests itself by 
the jump of the magnetization which is related to the appearance of abnormally large ±- 
contours. Without cutoffs this jump occurs for g ~ 1/N, and imposing the s(N) constraint 
would delay such a jump [3cS3j. It is easy to imagine what should be the critical order 
of the magnetic field g, at which those large contours should start to be favored in the 
s-restricted phase: for a db contour of the linear size s(N) one wins ~ s 2 g on the level of 
magnetization and loses ~ s on the level of surface energy. These two terms start to be 
comparable when sg ~ 1. Therefore no particular deviation from the classical behavior 
should be expected as far as gs(N) <C 1. We refer to [[Sc], where all these heuristic 
considerations have been made precise. 



3.2.2. Basic local estimate on the KlogN scale. Actually UlScj it is enough to con- 
sider only the basic K log iV-scale: 

Lemma 3.2.1 (fSc|). Assume that a sequence of numbers {&tv} satisfies 

b N \ogN 
lim — 5 — = 0. 

iV^oo N 2 

Then, on the basic scale s(N) = KlogN, the estimate 
[m a = {M A fx% + a N 

1 n 2 n 2 „», (3-2.3) 



holds uniformly in domains A G T>n and in a^ G M^n[0, b^], where X/3 is the susceptibility 
under the pure state 



3.2.3. Super-surface estimates in the restricted phases. Moderate deviations on 
the intermediate scales s(N) 3> logiV are, for the purposes of the theory, controlled by 
the following super-surface order estimate in the phase of small contours (c.f. Lemma 2.5.1 
in |S§) 
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Lemma 3.2.2. Let the large contour parameter s(N) S> logiV be fixed. There exists a 
constant c\ = c\((3) > 0, such that for all N > 0, A G V]y and all G 

M A = (M A f X :_ + a N ) < exp ( - A ). (3.2.4) 

The idea of the proof is simple: either an area of order aAr/2m* is exhausted by the 
X log AT large contours, which, in the n N _ -restricted phase, should have a surface ten- 
sion price with the exponent of the order a^/s{N), or if log N large contours cover an 
area much less than aN/2m*, which means that the remaining deficit of the magneti- 
zation should be compensated in the basic K log N restricted phase, where we can use 



Lemma 3.2.1 



3.3. Bulk Relaxation in Pure Phases 

The term relaxation is used here in the equilibrium setting in order to describe the 
approximation of local finite volume statistics by the infinite volume ones. We successively 
describe the relaxation properties of pure "— " states with non-positive and small positive 
magnetic fields and in the restricted phases of small contours. 

3.3.1. Non-positive magnetic fields h < 0. The crucial property of low temperature 
pure phases could be stated as follows: Let us say that the sites i and j are *-neighbors if 
\\i — j\\i = 1. Given a spin configuration a on { — 1, +1} Z , let us say that the sites i and j 
are +*-connected, if there exists a *-connected chain of sites i\, ...,i n , %\ = i and i n = j, 
such that a{ik) = 1 for every k = 1, ...,n. 

Theorem 3.3.1 ( |CCScj| ). For every (3 > (3 C there exists c\ = c\{0) > 0, such that 
uniformly in subsets A C Z 2 , i,j G A and in magnetic fields h < 0, 



A-,h ( * A j ) ^ e- clWll ^"°°. (3.3.1) 

Remark. Of course, since |i j\ is a non-decreasing event, the uniformity follows from 
the FKG ordering, once fl3.3.l|) is verified for the infinite volume zero-field measure /U_. 

Corollary 3.3.1 (Relaxation of local observables). Fix k £ 7L. Uniformly in A CI 1? , 
magnetic fields h < and local observables f with |supp(/)| = k, 

(f) A ,-,h ~ (f)-,h ^ C2 (A;)e- C3 ^ disto °( supp ( / ^ A ) (3.3.2) 

Furthermore, 

Corollary 3.3.2 (Relaxation and decay of semi-invariants). Fix n G Z. Uniformly in 
A C Z 2 , magnetic fields h<0 and sites i±, ...,i n G A, 



(a(i 1 )-..,,a(i n )f A _ h - (a(i 1 )-..,,a(i n )f_ h < Cl ( n ) e - c »03)« ltot -(^--WA) (3 3 3) 



and 



(ff(»i);..-;ff(»n)>x,-,h < c 6 {n)expi-c 7 {P) K - (3.3.4) 



Finally, 
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Corollary 3.3.3 (Asymptotic expansions). Fix n E Z. Uniformly in A C Z 2 and in 

i G A, 



k=l 



< c 8 (n)\h\ n+1 + C9 (n)e- CloW)dist °°( l > dA ) , 

(3.3.5) 

where is the k-th semi-invariant of the zero-field infinite volume measure y?_, 

Sk = ^2 (a(0);a{ii);...;a(i n )f_. 

ii,...,i^GZ 2 



Remark It is possible (and straightforward) to formulate ( [3.3.3| ), (|3.3.4| ) and ( 3.3.5| ) in 



the general case of n local observables fx, f n . □ 

3.3.2. Positive magnetic fields h > 0. Modifying "— " states by negative magnetic fields 
h < amounts to moving away from the phase transition region. Relaxation properties 
of yP A _ . with h > are radically different - uniformity is lost, and the size of the domain 
A starts to play a crucial role. Indeed, the unique infinite volume measure . = fx? 
stochastically dominates whatever small h > is. Thus, for large domains A, the 
configuration in the bulk is flipped under /j, A _ , into the "+" dominated state. It is easy 
to understand on the heuristic grounds what should be the order of the critical size of 
A for such a "flip" to occur: given h > 0, the surface energy of a ±-contour 7 is of the 
order I7I and it competes with the bulk gain inside the contour which, in its turn, is 
proportional to /iArea(7). The latter factor wins (loses), once the linear size of 7 is much 
larger (respectively much smaller) than \jh. Thus the sign of the dominant spin under 
I 1 A - h snom d depend on whether A can accommodate large enough contours, or, in other 
words, on how the linear size of A relates to 1/h. 

The important and remarkable fact is that exponential relaxation properties of finite 
volume "— " states are uniformly preserved for domains of the sub-critical size. 



Theorem 3.3.2 ( ffcSg| , |Sc|). There exists a constant a = a{(3) > such that for any 
h > fixed, 

P ( 4 4±U n \ < P -ci(/3)IK-j||oo 



A-,h ( * ) ^ e- ClWI|l - Jl|o °. (3.3.6) 

uniformly in domains A C Z such that any connected component of A has diameter 
bounded above by a/h. As a consequence exponential decay of semi-invariants ( [3.3.4 ) 



and the asymptotic expansion estimate ( 3.3.5| ) hold uniformly in such domains as well. 



3.3.3. Phases of small contours. Theorem |3.3.2| explains how the cutoff parameter 
s(N) upgrades the regular behavior of "—"-states with positive magnetic fields h: By the 
definition of the restricted phase \i A _ the diameter of any relevant microscopic domain is 
at most of the order s(N). 

Theorem 3.3.3 ( IgcSjj , |Sc|). There exists a constant a = a{(3) > such that for any 
h > and s satisfying hs ^ a(P), 

A S - h (i^j) < e-^W-'H- , (3.3.7) 
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uniformly in domains A C Z . 

Furthermore, the expectations in restricted phase are controlled as follows: for every k E Z, 

(/)!',-,, - </>inA s(/)> -, ft < ^K^, (3.3.8) 
uniformly in A C Z 2 and in Zoca/ functions f , | (supp(f)) | = fc, where we have used the 

following notation: A s (/) = {i : doo (i, supp(/)) ^ s}. Finally, the decay of the semi- 
invariants is controlled in the restricted phases as 

(a(i 1 );..-a(i n )f X ^ h \ < c 4 (n)exp j-c^ ^ 00 ^ 1 ' l A As j, (3.3.9) 

3.4. Calculus of Skeletons 

The renormalization analysis of large ± contours is performed on various cutoff scales 
s, the appropriate choice of s typically depending on the linear size N of the system 
s = s(N). We shall state coarse graining estimates uniformly in finite domains A C Z 2 
and in the cutoff scales s. 

3.4.1. Definition. A ± contour 7 is said to be s-large if diam oc (7) ^ s. Given a cutoff 
scale s E N and an s-large ± contour 7 we say that S = (u\, ...,u ra ) is an s-skeleton of 7, 
7 ~ S if 

1. All vertices of S lie on 7. 

2. s(N)/2 ^ — Uj+iHoo ^ 2s, V i = 1, ...,n, where we have identified n n+ i = ui. 

3. The Hausdorff distance c?h between 7 and the polygonal line Pol(S') through the 
vertices of S satisfies 

d H ( 7 ,Pol(S)) < s(N). 

Similarly, given the collection (71, 7„) of all sdarge contours of a configuration a E 

let us say that a collection & = (Si, S n ) of sdarge skeletons is compatible with a, a ~ ©, 

if 7i ~ Sj for all i = 1, n. 

Of course, a configuration <t E has, in general, many different compatible collec- 

tions of s-skeletons. Nonetheless, for each particular S the probability 

fi%_(6) i ^_(a: a~6) (3.4.1) 

is well defined. 



3.4.2. Energy estimate. As the renormalization scale s grows, the probabilities ( [3.4.1 ) 
start to admit a sharp characterization in terms of the energies Wp(&), 

n 

W P (&) t ^W^Pol^)), 
1 

for a collection © = (Si, S n ). Below we a give precise version of this crucial statement 
in terms of the upper and lower bounds on the corresponding probabilities. The first 



important renormalization energy estimates could be [Pf2] formulated as follows 
Lemma 3.4.1 (jPf|). On every skeleton scale s and independently of A C Z , 

S ) ^ exp{ -Wp(e) }. (3.4.2) 
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Furthermore, uniformly in A C Z ; r > and cutoff parameters s, 

ci log|Ak 



Ml_(W^(@)^r) <exp|-r(l 



(3.4.3) 



Energy estimate ( 3.4.2Q provides an upper bound on the probability of observing ± con- 
tours in the vicinity of a skeleton. Before going to a complementary lower bound let us 
dwell on the sample path structure of the contours which is hidden behind these renor- 
malization estimates. 



3.4.3. Calculus of skeletons. By definition a contour is a self-avoiding closed path of 
nearest neighbor bonds of Z 2 . For every set A C Z 2 the Ising measure yP A _ induces a 



weight function q p A * on the space of such self-avoiding polygons (see Subsection 1.2.2) 



Qa* (7) = - ( a e ^ : 7 is a =t contour of a) . 

In terms of these weights the probability of observing a certain skeleton S = {u%, ...,u n } 
could be written as 



E 

7-5 



J* 



(7). 



Each microscopic contour 7 compatible with S, 7 ~ S, splits into the union of disjoint 
open self-avoiding lattice paths 7& : — > Uk+i, k = l,...,n. The analysis of limit 
properties of fi A _ comprises two main steps which could be loosely described as follows: 
1) As the renormalization scale s grows, the statistical behavior of different pieces 7^ 
decouple under q A , , that is 



(3.4.4) 



7~S 



2) The k—th term (A; = 1, n) in the above product corresponds to a ± interface stretched 
in the direction of the vector Uk+i — Uk G M 2 , in other words 



9a* (7fc) 



(3.4.5) 



Thus, the skeleton calculus resembles a refined version of the sample path large devia- 
tion principle for genuinely two-dimensional random curves. At very low temperatures, a 



very precise local analysis of the phase separation line has been developed in | DKS |,|DS| 
using the method of cluster expansions. Our approach here pertains to the whole of the 
phase transition region ft > f3 c , but is strongly linked to the very specific self-duality prop- 



erties of the two-dimensional nearest neighbor Ising model. We refer to Subsection 1.2.2 



and, eventually, to [PV2, PV3] for comprehensive description and study of the relevant 
properties of the duality transformation. The output of these techniques could be recorded 
in the following form 



Lemma 3.4.2 (Probabilistic Structure of the Phase Separation Line | PV2|| ). Given any A C 
1? and any two compatible self-avoiding paths Ai and A2, 



4*(AiUA 2 ) ^<4(Ai)<&(A 2 ). 



(3.4.6) 
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Furthermore, 



-ci(/3)|A 2 | 



g£(AiL)A 2 ) 
q% (Ai) 



-C2 03)|Aa| 



(3.4.7) 



On the other hand, given any AQX 2 and any three points u,v,w £ A*, the q 1 ^* weight of 



the paths going from u to v through w is bounded above as [ P V2 ] 



E w < I E £ (AO E £ (a 2 ) ) • 

Ai:u— »ui / \A2:ui— *v 



(3.4.8) 



A:u— >?; 
uieA 



Finally, the weights q A * are non-increasing in A, and are related to the dual connectivities 



as 



E d*(A) = m«m«)>s; 



(3.4.9) 



A: u — >v 



Relation ( 3. 4. 9] ) is the link to the surface tension: first of all the impact of a particular 
set A exponentially diminishes with the distance to dA EL 

(a(u)a(v))f -exp{-c 2 ((3)d({u,v},dA)} < {v(u)a(v))%.j < (a(u)a(t;))f . 

(3.4.10) 

uniformly in ^4* C Z 2 and any u, v G A*. Moreover the following Ornstein-Zernike type 
correction formula jAj holds uniformly in u, v G Z 2 : 

exp{-T0 (« - w) - c 3 (/?)log ||«- «!!«,} < {o-{u)a{v)Y f ^ exp {-Tp (u - v)} , 

(3.4.11) 



3.4.4. Skeleton lower bound. The energy estimate (3.4.2) is an immediate consequence 
of the (iterated) sub-multiplicative property (|3.4.8|) , the representation formula ( 3.4.9j ) and 
the right-most inequalities in ( |3.4.10| ) and ( |3.4.11|). In order to prove a lower bound one 
essentially needs to reverse the inequality in (|3,4.8D. An indirect way to do so is to use 



the FK representation (see [ |5cS2 | and |ISd] ). We shall briefly present here a more direct 
approach which has been developed in [M and |PV2 ] . Qualitatively it gives the same order 
of corrections as the FK one, but has a clear advantage of being explicitly related to the 
statistics of the microscopic phase boundaries at different length scales. The basic idea 
is that the phase separation line has rather strong mixing properties, in particular paths 
Ai and A2 on the right hand side of ( |3.4.8 ) should interfere, in the case of (u,v,w) being 
in a general position, only in a vicinity of w. Thus, at a price of lower order corrections 
(as we shall see these corrections are logarithmic with the skeleton scale s) the inequality 
( 3.4.8 ) could be reversed using the super-multiplicativity property (|3.4.6| ). The notion of 



"general position" simply means that u, w and v do not form too small an angle and live 
on the same length scale, and it is quantified by the following 

Definition. Given a skeleton scale s £ N and a number e > 0, let us say that that a triple 
(u, w, v) of Z 2 -lattice points is (s, ^-compatible, if 



— ^ min {\\w — u\\ 



\v-w\\ x } ^ max{\\w - u\\ x ,\\v - wW^} < 2s, 



whereas cos (w — u,v — w) ^ — 1 + e. 

We shall state the lower bound in terms of the limiting weights (•) = lim^* 
(which exist by Lemma 3.4.2). 



□ 
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Lemma 3.4.3. Fix e > 0. Then there exists a scale s = s(e), such that 

/*(A) > exp{-(Tp(w-u) + T (v-w)) -ci 09) logs}, (3.4.12) 



X:u— >v 



uniformly in all skeleton scales s s(e) and in a/l (s,e)- compatible triples (u,w,v). 



We sketch the proof of this lemma in Appendix B. Iterating ( 3.4.12 ) we arrive to the 
following lower bound on the probability of observing a certain regular skeleton: 
Definition. A skeleton S = (ui, ■■■,u n ) is said to be (s, e)-regular, if any triple (itj_i, Uj, Ui+i) 
of successive points of S is (s, e)-compatible, and the distance between any two non- 
neighboring intervals [uj,Uj + i] and [v,j,Uj+i] exceeds es. □ 

Lemma 3.4.4. For every e > 0, there exists a number C2 = 02(e) < 00, such that uni- 
formly in the skeleton scales s and in all (s , e) -regular skeletons S, 

/%- (3 a ± contour 7 : de(7, Pol(5)) < K(j3) \/slog s) 

> exp{-^(Pol(5)) - c 2 (e)#(S)logs} {3A U) 



> exp|-^(Pol(5))(l - c 3 ( £ ,/3) l0gS 



where #(S) denotes the number of vertices in S, and the last inequality follows from 

In fact we need lower bounds only for a very specific set of s-skeletons, namely on those 
approximating the Wulff shape tC aN /2m*- These skeletons always satisfy the conditions of 
the above theorem. An academic attempt to prove a lower bound for all possible shapes 
will lead to annoying, though solvable, technicalities, but will fail to contribute much to 
the microscopic theory of phase separation, as we see it. 



3.5. Structure of The Proof 

In order to give a probabilistic characterization of the microscopic canonical state 
H N _ ( ■ |Mtv = —m*N 2 + cin) one first derives a sharpest possible lower bound on the 

probability fi N _ (Mjy = —m*N 2 + aN), and then rules out those geometric events (in 
terms of skeletons, but with an eventual translation to the language of microscopic spin 
configurations), which happen to qualify as improbable when compared with this lower 
bound. 



3.5.1. Lower bound. The best lower bound comes as an outcome of the optimal com- 
bination of the basic local limit Lemma 3.2.1 and the skeleton lower bound ( |3.4.13| ). We 
choose a skeleton approximation of the corresponding Wulff shape /C ajv /2m*i an d using 
local limit estimates steer the magnetization towards the desirable value —m*N 2 + a^. 
Optimality reflects the choice of the best possible skeleton scale: Notice that the estimate 
(3.4.13) becomes sharper with the growth of the cutoff parameter s(N). On the other 
hand, the area of the microscopic phase region is controlled, with respect to the area in- 
side Pol(S) 



ajv/2m*, up to a iVy s(N) log s(N) correction (see Appendix B or pSdJ ), 
which, of course, makes the local limit step more expensive for large values of s(N). It 
happens that the bounds are balanced on the skeleton scale s(N) ~ ^/on- 
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Theorem 3.5.1 ([ [[Sc| ). Uniformly in € that is for all aj^ 0, such that 

-m*N 2 + a N <E Range(MTv), 

Mat,- (Mjv = -m*N 2 + a N ) > exp j-y^y^ {dKL x ) - Cl (/3) ^ log ivj . (3.5.1) 

3.5.2. Upper bounds. First of all, one derives an upper bound on the shift of the mag- 
netization. On any skeleton scale, 

Mat,- (Mjy = -m*N 2 + a N ) < ^ _ (Mjv = -m*iV 2 + ; 6) . (3.5.2) 

e 

Due to the intrinsic entropy cancelation under the skeleton coarse graining, and in view 
of the lower bound ( |3.5.1 ) and the energy estimate ( |3.4.2j ) one could, for example, shoot 



for the maximal term in the above sum. If the phase volume (sec [ DKSfl for the precise 



definition ) of 6 is much less than aj\[/2m*, then the deficit of the magnetization should be 



compensated in the phase of s(iV)-small contours, which, by Lemma 3.2.2 exerts a super- 
surface price in the exponent. On the other hand, if the phase volume of © is close to 
a 7v /2m* , then by the isoperimetric inequality and by the energy estimate 

(13.4.2| ), the best 

possible price one should be prepared to pay is already close to expj — Wp (t£ aN /2m*)}- 
Again the resulting estimate is subject to an optimization via a careful choice of the 
skeleton scale s(N). 

Theorem 3.5.2 (pep. Uniformly in a N ~ N 2 , 
A,- 



(M N = -m*N 2 + a N ) < exp l-^LyVp (did) + ci(/?)^logivJ . (3.5.3) 



A more delicate study JDKS| ] , [ ISc ] of the typical sample properties of the microscopic 
configuration a under pL N _ ( • \Mn = —m*N 2 + a at) is again based on the analysis of 
(]3 . 5 . 2[) . At this point the stability Bonnesen-type estimates (see Subsection 1.3 of the 
Introduction) for the Wulff variational problem become important - they enable to quantify 
the conclusion that only those collections S, which are close to the shifts of the Wulff 
shape /C ajv /2m*> have a chance to survive a comparison with the lower bound ( 3.5. 1| ). A 



step further, involving local limit estimates of Lemma 3.2.1, is to conclude that all these 
collections actually contain exactly one large skeleton, which corresponds to the unique 
large contour as asserted by the DKS theorem. 

3.6. Open Problems 

There are still important open problems even in the nearest neighbor Ising case. No- 
tably, one knows how to control precise fluctuations of the phase separation line only at 



very low temperatures, that is using the method of cluster expansions [DH]. This is a 
serious gap in the theory, since large scale statistics of microscopic phase boundaries are 
ultimately responsible for exact (up to zero order terms) expansions of canonical partition 
functions Q. So far qualitative probabilistic results have been obtained either for very 
low temperature models [Q, or in the simplified setting of self- avoiding polygons p3[| , [HI] 
or Bernoulli bond percolation [CI|. Another interesting and apparently important prob- 
lem is to understand sample path properties of spin configurations in a situation when 
a canonical constraint is imposed in the restricted phase. Apart from giving rise to a 
potentially fascinating probabilistic structure, this question is closely related to the issue 
of the dynamical spinodal decomposition. 
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There is absolutely no matching probabilistic study of the phase separation in multi- 
phase two-dimensional models, for example (/-states Potts models. Some results in this 
direction are reported in Q, but this issue is almost entirely open even in the context of 
the Li-theory. In particular, the corresponding phenomena is still not worked out on the 
level of macroscopic variational problems, see, however [ABFH], [MoS| and the references 
therein. 

The key issue, however, which we feel is largely misunderstood is that at moderately 
low temperatures the DKS theory of two-dimensional phase segregation, say in the gen- 
eral context of finite range ferromagnetic models with pair interactions is far from being 
complete. What currently exists is an example of how these ideas could be implemented 
in the nearest neighbor case. At least from the mathematical point of view, the nearest 
neighbor case is a degenerate one, in a sense that it enables a reduction to pure bound- 
ary conditions over decoupled microscopic regions even at temperatures only moderately 
below critical. This should not be the case for more general range of interactions. In 
this respect the assertion that low temperature expansions should go through for general 
interactions much along the same lines as they do for the nearest neighbor model, seems 
to be rather irrelevant - the real issue is not to kill mixed boundary conditions, but to 
understand how they should be incorporated into the DKS theory. 
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Part 4. Boundary effects 

In the previous parts, we explained how the thermodynamical variational problem de- 
scribing the macroscopic geometry of coexisting phases can be derived in various lattice 
models of statistical physics. To simplify the analysis, we restricted our attention to peri- 
odic boundary conditions or to systems contained in a Wulff-shaped box, avoiding thus a 
discussion of the effect of a confining geometry on the behavior of the system. In this part, 
we would like to explain what happens when we take such effects into account. Boundary 
conditions play a particularly important role in the kind of problems presented in this 
review, since they concern the asymptotic behavior of large but finite systems and there- 
fore the boundary cannot be simply "sent to infinity" as usually done. We will see that 
taking care of boundary effects not only provides a complete description of the geometry of 
these constrained systems thus allowing a rigorous description of the interaction between 
an equilibrium crystal and a substrate, but also allows to study the effect of so-called 
boundary phase transitions. 

For simplicity, we only discuss the case of the Ising model with nearest neighbors inter- 
action. 



4.1. Wall free energy 

The vessel containing the system has not only the property of confining it, but can also 
act in an asymmetric way on the various phases inside, favoring some of them; indeed 
this is what happens typically in real systems. In fact, this is precisely the reason one 
introduces boundary conditions in the first place: To impose the equilibrium phase the 
system realizes. It appears to be convenient to have a parameter allowing a fine-tuning 
of the asymmetry, interpolating between pure + or — boundary conditions. Let us now 
describe how this is done. 

Let £ = {i G Z d : i(d) = 0} and h d = {i G Z d : i(d) ^ 0}. The vessel of our system is 
the box 

B NM = {i G h d : -N i(n) ^ N, n = 1, . . . , d - 1, ^ i{d) < M} , 

and the wall is Sat = Djv,M n E. 

Let n G M; we consider the following Hamiltonian, 

H o N>M ( a ) = - Yl ai<y i - * Gi ■ 

(i,j)Ch d iGSjv 

Let a G {—1, l} Ld ; the Gibbs measure in Btv,m with boundary condition a is the following 
probability measure on { — 1, l} L<i [], 

n,m,& 1q otherwise. 

We'll usually use the short-hand notations ^ v Mw , Z^^-, .... As usual, we write + 
for a = 1 and — for a = — 1. We therefore distinguish one of the sides of the box 

9 Note that we could equivalently consider fj^' 71 _ as a probability measure on { — 1, 1} Z by extending 
the b.c. a by <fj = 1 for all i 6 Z d \ h d ; it is then possible to replace the boundary magnetic field n by a 
coupling constant: ^X^eEjv ff i = V S(j j)- ies N j£L d a i a i- This will be used when dealing with negative 



boundary field, see Subsection 4.4.1 
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Bn^m, ^n, which we call the "wall". Notice that instead of usual boundary conditions, 
a boundary magnetic field r\ is acting on S^r; since setting 77 = 1 produces + b.c. on 
the wall, while setting 77 = —1 results in — b.c, this provides the promised interpolation 
parameter. Of course, we could also consider more complicated situations, where (possibly 
inhomogeneous) boundary magnetic fields act on the whole boundary of the box. However, 
for simplicity, we restrict our attention to this particular case, which will turn out to be 
general enough that the basic phenomena induced by the use of boundary fields can already 
be analyzed. 

To quantify the preference of the wall toward one of the phases, it is convenient to 
introduce a new thermodynamic quantity, the wall free energy, 

1 7?" 

r bd (A 77) ± lim -— log -»± . (4.1.1) 
IEjvI Z^ M> _ 

The existence of this quantity, and the remarkable fact that the two limits can be taken 
in any order, has been established in [FP1|; the proof relies on the simple identity 

T hd {M = Hm P f d V- (4.1-2) 

We'll return to this formula in the next section. The heuristics behind the definition of 
T bd{P,v) i s that the free energy F^ V M + ,_^ = — log Z^^- of the + (— ) phase can be 
decomposed in the following way: 

F N%,+ = MP) \®>n,m\ + ftiP) \9®>n,m \ Siv| + / W + (M |£iv| + o(\dB N ,M\,\Z N \) , 

F^ M> _ = W) \B N>M \ + frifi) \d® N ,M \ ^n\ + f~(P,v) |Sjv| + o(\dB NM \,\E N \) , 

where 

/ b (/3) = - lim iD^Ml-MogZ^-, 

/+(/?) = - lim laO^I -1 (log Zfu + - fM\®N,M\) > 

i\.,Ad — too ' ' 

/ W + (M ±- N KmjZ N r 1 ( lo g Z K, M ,* " W)\®n,m\ ~ /+(/?) I \ E*|) , 

(and similarly for f~((3) and /~ (/3, 77)). As the notations suggest, f b ((3) is independent 
of 77 and a, / s + (/3) is independent of 77 and by symmetry / a + (/3) = / 8 - (/3). Therefore, we 
see that r bd (/3, 77) = lim^oc t^-t (Ffi v M _ - Ffr' v M ,) = f~(/3, 77) - f+(/3, 77) is nothing else 

than the leading order term of the difference in free energy between the two phases in the 
presence of the wall. 

The ultimate justification of ( |4.1.1| ) however is that this quantity plays exactly the role 
of its thermodynamical analogue in the variational problem describing the macroscopic 
geometry of phase coexistence, see Theorems 4.3.2 and [4.3.3] below. 

The following Theorem states basic properties of r bd (/3, 77); since r bd (/3, 77) is obviously 
odd in 77, we just state them for 77 ^ (also r bd (/3, 0) =0). 



Theorem 4.1.1. | FP2 | Let = Tg(e<i) and suppose 77 ^ 0. Then 

• r bd (/3, 77) is a non-negative, increasing function of and 77, concave in 77; moreover, 
if 77 > 0, 

(/3, 77) >0^/?>/3 c . 
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For all (3 and rj, T hd (f3,rj) ^ t%. 

For all (3 > (3 cl there exists 1 ^ ?7 w (/3) > such that 

% rj) <T*p^r]< Vw ((3). 



In the case of the 2D Ising model, ^((3) can be computed explicitly, see [ Abl , MW ] 
and Fig. |n|. 

The following terminology is standard^ when n f?„(/3), we say that the system is 
in the complete drying regime; when \rj\ < rj w (f3), it is in the partial wetting regime; and 
when j) ^ — r/ w (/3), it is in the complete wetting regime. The reason for this terminology 
should become clear later. 

4.2. Surface phase transition 

In this section, we will see that the boundary magnetic field can trigger surface phase 
transitions: The behavior of the system in the vicinity of the wall depends dramatically 
on | r/\ being greater or smaller than r/ w (/3). A more detailed discussion of these issues can 



be found in PV1 ], 



The state of the system in the middle of a big box Btv,m is entirely determined by 
the boundary conditions, and is independent of the value of the boundary field, so that 
the usual (infinite volume) Gibbs state simply doesn't provide any information on the 
behavior of the system close to the wall. To analyze the behavior of the system "in the 
vicinity" of the wall, it is therefore useful to introduce the notion of surface Gibbs states; 
these differ from the Gibbs states usually considered in these models by the fact that 
one does not work with a sequence of boxes converging to but instead converging 
only to the half-space h d . More precisely, the surface Gibbs states are the weak limits of 
the measures /^Ma when N, M —* oo (observe that B/v,M S L d ). Two of them are of 
particular importance for our discussion, fJ^a + and obtained respectively by taking 

weak limits of the measures with + and — boundary conditions. It is not difficult to show 



[FP1] that these two measures exist, are extremal, and are invariant under translations 
parallel to the wall; moreover, there is uniqueness of the surface Gibbs state if and only if 

There is a close relation between r bd (/3, rj) and the behavior of the system near the 
wall; this can be most easily seen from the following identity, consequence of ( 4.1.2; ) and 
symmetry [ FP1[| , 



r bd (M = / (°i>)& + di/ = / (Wofj' + ~ (vofj'j dr/. (4.2.1) 



Using (4.2.1), it is possible to prove the following Theorem showing that a surface phase 



transition occurs at rj = r] w ((3); this is the so-called wetting transition. 



Theorem 4.2.1. [FP2] There is a unique surface Gibbs state if and only if |ry| ^ rj w (f3). 



Let us briefly discuss the heuristics behind this result. The + and — boundary conditions 
fix the phase present in the bulk (i.e. in the middle of a big box Ojvm)- However, Theorem 



4.2.1 shows that when rj ^ t? w (/3), the surface Gibbs state is unique, and therefore the state 



of the system near the wall is independent of the boundary conditions, i.e. of the phase 



10 This terminology only makes sense once we have chosen one of the equilibrium phase as reference; 
here it is the — phase. 
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present in the bulk. The mechanism responsible for this is the following. Suppose that 
77 < and consider +-boundary conditions; then it is natural to regard the boundary 
field as a negative b.c, and therefore to introduce an open contour with boundary dY,^ 
separating the — phase favored by the wall from the + phase present in the bulk (see 



Section 4.4 for more details). As long as r\ > —1, there is a competition between two 
effects: On the one hand it is energetically favorable for the open contour to follow the 
wall, on the other hand this would lead to a loss in entropy, since there is less room 
for fluctuations. When r\ ^ — ry w (/3), the entropy wins: The contour is repelled away 
from the wall, at a distance diverging with the size of the box; this is the phenomenon 
of entropic repulsion. The surface Gibbs state then describes the behavior of the system 
below this surface, i.e. a mesoscopic film of — phase along the bottom wall. The fact 
that the contour is sent away from the wall explains why we recover the surface tension, 
T bd(P,v) = T p- When rj > — 77 w (/3) energy wins, and this modifies completely the behavior 
of the microscopic surface: it sticks to the wall, making only small excursions away from it; 
in this case, the phase in the bulk can reach the wall and the surface Gibbs state depends 
on the choice of boundary conditions. 

Part of these heuristics can be made quite precise in the 2D case. Consider + boundary 
conditions. When > rj > — r/ w (/?), one can prove that the probability that a connected 
piece / of the wall is not touched by the open contour is bounded above by K exp[— (rj — 



T bd(P,v)) l-^lh showing that the phase separation line really sticks to the wall [ PV2| ], The 



informations available when 77 ^ — r/ w (/3) are much less precise; the magnetization profile 



computed in [Abl] shows that there is a film of width of order y/ N along the wall. A 
related, much more precise result, which holds at sufficiently low temperature and for 
77 = —1 is that the phase separation line, once suitably rescaled, converges weakly to the 
Brownian excursion Q ; this should be true for any 77 ^ — t/ w (/3) . 

In higher dimensions, much less is known. When 77 > — r/ w (/3), one can show that the 
probability that the open contour touches the middle of the wall is bounded away from 



uniformly in the size of the box | FP2 |. When 77 ^ — t/„(/3), very little is known, except 
in the simpler case of SOS models. Also, if it is known in dimension 2 that r/ w (/3) < 1 
(since the exact expression for r? w (/3) has been computed | Ablfl ), this is an open problem 
in higher dimensions. 



Theorem 4.2.1 gives a first explanation of the terminology introduced above: when the 
system is in the complete drying regime, the equilibrium phase along the wall is the + 
phase, whatever the phase in the bulk is; when there is complete wetting, it is the — phase; 
only in the regime of partial wetting can both phases be present near the wall. The fact 
that the phase transition is determined by r/ w (/3) (i.e. the characterization of the partial 
wetting regime by = \r hd ((3, r/)\) is known as Cahn's criterion. 

4.3. Derivation of the Winterbottom construction 

In this section, we show how Winterbottom construction, describing the equilibrium 
shape of a crystal in the presence of an attractive substrate, can be recovered from a mi- 
croscopic theory. To do this, we consider the measure ^ v rN , , for some r£l, conditioned 
with some canonical constraint (exact or approximate, see below). Of course, the situation 
here is more complicated than the one described in the introduction, since instead of an 
infinite wall, the system is contained in a finite vessel. This, of course, makes the problem 
more difficult: When the solution of the Winterbottom variational problem does not fit 

inside the box Df = {1 £ l d : \x(n)\ ^ 1, n = 1, . . . , d — 1, ^ x(d) ^ r}, the solution of 
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Figure 11. The case of the 2D Ising model. Left: The phase diagram; the region of 
non- uniqueness of the surface Gibbs state is shaded. In the other region, there is a single 
surface Gibbs state. Right: A sequence of equilibrium shapes. 



the constrained problem will differ from Winterbottom shape. In fact, the general solution 
of the constrained problem is not known. In the way we state them below, the derivation 
of this variational problem from statistical mechanics still applies in the case when the 
solution is not known. 

Before stating the main Theorems of this Part, we briefly describe how the wetting 
transition manifests itself in the macroscopic geometry of phase separation. To do this, 
let f5 > f3 c be fixed, and choose a value m for the canonical constraint so that the corre- 
sponding Wulff shape is small enough to be placed inside the box OJ°. If 77 ^ ?? w (/3), then 
T bd (/3, 77) = Tg, and the typical configurations will consist of a macroscopic droplet of — 
phase, with Wulff shape, immersed in a background of + phase; in particular, the shape 
of the droplet is independent of the value of the boundary field (Fig. [IT] a). This behavior 
persists up to the value 77 = r] w ((3). Notice that as soon as 77 < 1, it becomes energetically 
more favorable for the droplet to touch the wall. In dimension 2, however, since 7y w (/?) < 1, 
the droplet stays away from the wall, because entropy loss is not compensated by energy 
gain until 77 reaches the value 77 w (/3). It is an interesting open problem to decide whether 
WwiP) = 1 for d > 2. When 77 < 7? w (/3), the typical configurations consist of a macroscopic 
droplet, with Winterbottom shape, tied to the wall. The shape of the droplet now depends 
on the value of 77, and decreasing the boundary field amounts to letting the droplet spread 
more and more (Fig. b-e). For some value 77, the droplet covers for the first time the 
entire wall (Fig. |ll| e). From this point on, the shape of the droplet is left unchanged when 
77 is decreased (Fig. |ll| f; the dashed line represent part of a possible "true" equilibrium 
shape for the unconstrained problem). 

From this discussion, we see that the wetting transition at rj w ([3) has a macroscopic 
manifestation in the canonical ensemble. Because of the confined geometry, however, the 
second transition, at 77 = — t/ w (/3) cannot be seen. To be able to detect it, one has to 
consider mesoscopic droplets (in the form of large moderate deviations, see the remark 
after Theorem 4.3. 2| ) . 



This also explains pretty well the terminology introduced previously: In the complete 
drying regime, the droplet stays away from the wall, and so the wall is completely dry 
w.r.t. the — phase; in the partial wetting regime, the droplet touches the wall, and both 
the + and — phase are in contact with it (provided 77 < 77). The complete wetting regime 
cannot be distinguished from the partial wetting regime in this setting, but see the remark 



after Theorem 4.3.2 for a discussion of this issue. 
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4.3.1. 2D Ising model. Let r G 

configurations under the measure 



r'N.rN, 



where m G 



-771 



writing simply /-tjy?. 



The aim of this subsection is to describe the typical 

,[-\M N = m\B N>rN \), 

t (jj; we will simplify the notations further by 
As in Part ||[ it is possible to obtain precise 



m*) and M N = Z l&N , 
(r being kept fixed). 

asymptotics for the large deviations, in the form of the following generalization of the first 
part of Theorem 3.1.1 . Let W^ r; (m) be the infimum of the functional W/^ on subsets of 

ff,2| 



i 2 . with volume 



2m* 



Theorem 4.3.1. Let the inverse temperature (3 > (3 C and the boundary magnetic field 
77 G M be fixed; let the sequence {a^}; —m * \Bn,tn\ + a^v G Range(MAr), be such that the 
limit 



lim • 

N^oo 



aN 



»N,rN\ 



G (0,2m* 03)) 



exists. Then, 



log if?? 1 . 



{M N 
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T>N,r N \ ~ a N ) = -Wl v (1 + OiN- 1 ' 2 log N)). 



A version of this Theorem, in an approximate canonical ensemble (as in ( |4.3.1| )), has 
been proven in [PV2]; this stronger version can be obtained by combining the techniques 
of fvj and o f feci see Section fO. 

In Theorem |4.3.l| , we have made no statement about the asymptotic description of the 
typical configurations under the conditioned measure. The reason is the following: These 
strong concentration results require the knowledge of stability properties of the variational 
problem in the form, for example, of Bonnesen inequality. However, in the present case, one 
does not always have that much information about the variational problem; in fact, even 
its solution is not always known. This prevents us from translating the energy estimates 
on the skeletons (see ( |4.4.8[) , ( [4.4.10 ) and ( [4.4.11| )) into strong concentration properties 
of the microscopic contours. Of course, in the situations when such stability properties 
are known ( fKP| contains a simple derivation of such a result for many situations), it is 
possible to obtain statements of the same kind as those of Part ||. 

This illustrates the fact that although the probabilistic theory in the 2D case is com- 
plete, in the sense that all the relevant information on the microscopic scale is available, 
the sharpness of the statements one can make on the macroscopic scale still depends 
on macroscopic stability properties, which are logically separated from the probabilistic 
aspect of the analysis. 

However, even without information about the stability properties of the variational 
problem, it is still possible to derive weak concentration properties, in a Li setting close 
to the one of Part ||. We present such a result in the way it is stated in [ PV2| . In this paper, 
an approximate canonical ensemble was considered, i.e. the measure was ■ | -4(77i; c)), 

where 



A(m; c) = jer : 



»N,rN\ 



771 



(4.3.1) 



with —771* < 77i < ?7i*, and c is some real number not too large (see Theorem 4.3.2| below). 
We are going to prove that the phases concentrate near macroscopic droplets which belong 
to the set Vim) 



V{m) = \ V C W r : \V\ 



771 



771 



2?77* 



»} 
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Recall that to each V G T>(m), we associate the function ly = lyc — ly. 

To state this phase segregation Theorem, we use the mesoscopic notation introduced in 
Part ^. Recall that N = 2 n . For any a < 1, we define a magnetization profile Ai\ an \(o~,x) 

at the 2t an l-scale which is piecewise constant on boxes M n _^ an ^(x) with x E T^_[ an ]> 

M [an] (a,x)=2- d ^ ( 43 - 2 ) 

We get 

Theorem 4.3.2. |PV2(| Let f3 > f3 c , r) € R, -m* < m < m* and 1/4 > c > 0. T/ierc 
i/iere exist a function 5(N) such that liniAr^oo <5(iV) = 0, a rea/ number k > and a 
coarse- graining parameter 1 > a > suc/i i/ia£ /or iV /arge enough 

|J V(Iy,5(JV))| A(m;c)) ^ 1 - exp{-0(iV K )} • 

\/eD(m) 



Remark: In this case, it should also be possible to study the whole range of moderate 
deviations, combining the techniques of [ISc] and [PV2], although this has not been done 
explicitly. We briefly describe the results obtained for large deviations sufficiently close to 
volume order jv|]. 

As long as r\ > — rj w (j3), the results are similar to those obtained in the setting of Part ||: 
The measure concentrates on configurations containing a single large droplet of — phase, 
with Wulff or Winterbottom shape depending on rj; in particular, the order of the large 
moderate deviations is still exp{— 0(y / a]v)}- There should not be any problem to extend 
this to the whole large deviations regime (a^ ^> AT 4 / 3 ). 

More interesting is the case rj ^ — ry w (/3). For those values of the boundary field, the 
system is in the complete wetting regime (r bd (/3, 77) = —Tp), and the solution of the un- 
constrained variational problem is degenerate. The solution of the constrained variational 
problem in is however still well-defined for every iV; it is obtained by extracting the cap 
of a Wulff shape and rescaling it so that the basis of the cap completely covers the wall and 
the rescaled cap has the required volume. When N goes to infinity, this droplet spreads 
out to become a thin film in the limit (covering the entire wall, hence the terminology 
complete wetting) , and the corresponding value of the surface free energy functional goes 
to zero. As a result of this, the scale of the large moderate deviations is not the same 
as when rj < r/ w (/3); indeed the leading term of the asymptotics can again be computed 
explicitly, and is found to be of order exp{— 0{{a^) 2 _/V~ 3 )}. In particular, we see that the 
large moderate deviations cannot extend up to ajq ~ 

AT4/3 

since (ajv) 2 N 3 is of order 1 

already when ~ iV 3 / 2 . This should not be surprising since, in the complete wetting 
regime, the volume under the microscopic contour is expected to have typical fluctuations 
of order N 3 / 2 (this can be shown when rj = — 1 and (3 is very large using the convergence 
to Brownian excursion stated in 0). Therefore, typical fluctuations of magnetization 
in the complete wetting regime are not governed by bulk fluctuations anymore, but by 
fluctuations of the microscopic phase separation line. To prove that this behavior is valid 
up to on ~ iV 3 / 2 might be a non-trivial task. □ 
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4.3.2. Ising model in D ^ 3. Let r£R and let T>(m) be the set of macroscopic droplets 
at equilibrium in B^f, 

V{m) = {ycif: \V\ = , Wp, v {dV) = W£„(ro)} . 

The rest of the notations were introduced in Part ||[ The main result is the following 



Theorem 4.3.3. pTV| For any (5 in 53 p; any rj G R, any m in (— m*,m*), the following 
holds: For any 5 > 0, there is ko = ko(6) such that for v < k 



lim min 

N^co k ^ k sg vn Jv > 4 



U V(lv,fl 

\ VeV(m) 



Mm ^ m |D7VrAf| 



VeD(m) 

4.4. The tools 



In this Section, we explain how the procedures described in Parts Q and |3| have to be 
modified to take into account the effect of the boundary. 

4.4.1. 2D Ising model. We describe the main modifications one needs to apply to the 
proofs of Part |3] in order to get the results stated in Theorems 4.3.1 and 14.3. 2j . We split 



this Subsection into two parts, one dealing with the lower bound on n^_(A(m; c)) or 
I1^P_{Mn = — m*\I]>N,rN\ + (in), the other one with the upper bound. 

The lower bound. The constrained variational problem is more difficult than the usual one. 
In fact, as noted above, the solution (and a fortiori its stability) is not known in general, 
although it is in many cases. This prevents us from proceeding as in Part ||, where the lower 
bound follows from summing over large contours fluctuating around the Wulff shape. It 
would then appear necessary to make the same kind of proof, but for any configurations of 
droplets surrounding the right volume (all potential solutions to the variational problem). 
This, however, would be tricky; indeed, since we want our results to hold for large, but 
finite boxes, it is compulsory to obtain estimates uniform over the droplet in the chosen set! 
Fortunately, properties of the surface tension and wall free energy allow us to restrict our 
analysis to a small class of well-behaved droplets: The solution of the variational problem 
is necessarily taken on a single convex droplet. This is a consequence of the convexity of 
Tp (use Jensen inequality) and the fact that r bd (/3, rj) ^ r|, which imply that replacing a 
droplet by its convex hull cannot increase the surface free energy; rescaling the resulting 
droplet decreases the energy even more. It is thus enough to prove the following 



Proposition 4.4.1. [ |PV2| Let [3 > (3 C and T] G R. There exists Nq = No(j3,rj,m,c,r) and 



a constant C such that, for any simple closed rectifiable curve C which is the boundary of 
a convex body of volume |Bj?|(m*(/3) + m)/2m*(/3) contained in and for all N ^ iVo, 

^ v _(A(m;c)) ^ exp{-W A „(C) N-fiCN 1 / 2 log N} . 
A completely analogous statement holds in the case of the exact canonical ensemble. 



The proof of Proposition [4.4. 1| is similar to the proof of Theorem |3.5.1| . We sketch now 
the main changes needed to deal with the boundary conditions. The case rj ^ requires 
a slightly more complicated proof than the case rj > so we first consider the latter. 

First case: n > 

As in the usual case, we want to approximate C with some polygonal curve with vertices 
on the dual lattice, and then sum over all contours going through the latter; this would 
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Figure 12. When Tbd(/3, »?) < Tp, the open contour connecting two sites close enough 
to the wall might not stay inside an elliptical set as in the bulk (dashed contour), but 
instead might get pinned by the wall (full contour). In such a case, the exponential 
decay-rate is in general not given by Tp or r bd (/3, rj) . 

allow us to extract, for each piece of the contour, the surface tension of the corresponding 
part of the polygonal line. Here, however, we want to be able to extract the wall free 
energy when the curve C follows the wall. There are some complications related to this: If 
two vertices are close to the wall, but don't belong to itPl, the sum over the corresponding 
piece of contour might not yield simply or r bd (/3, if), but some complicated mixture, 
since typical such contours might first go down to the wall, then follow it on some length, 
and only then go up to the other vertex, see Fig. Il2| this kind of behavior has been studied 



in details in [PV3]. It turns out that it is possible to construct a polygonal approximation 
to the curve C whose surface tension is not too large in comparison with that of C, while 
removing these possible pathologies. 

The idea is the following. Let 5n = N~ 1 / 2 log N, and set 

B>l(N) = {i 6 §J : min \\y - x\\ x > 5 N } . 

Let V be the convex body with boundary C and set Cn = d(V DB^(iV)) We first construct 
a polygonal approximation for each of the components of Cn n B^(iV) with segments of 
length 5n (apart from at most 8 of them which may be shorter). Set [x,y] = {z S Cn ■ 
z{2) = 5n}- If [x, y] 7^ 0, we connect the two corresponding pieces of polygonal lines by a 
broken line from x to (x(l),0), then to (y(l),0), and finally to y; we divide the segment 
between (x(l),0) and (y(l),0) into segments of length 5at/2 (except possibly for the last 
one which can be shorter). We repeat this construction for the three other sides of the 
box. The resulting closed polygonal line is denoted by Vn (see Fig. |l3|). Notice that by 
construction there exists an absolute constant C such that 

|vol(C) -voYP N \ ^C\%l\8 N . 

We then rescale the polygonal line Vn by a factor N and if necessary move slightly the 
rescaled vertices so that they belong to the dual lattice; the rescaled polygons is denoted 
by Vn- We then define a class (5 of closed contours going through the vertices of Vn 



11 Consider, for example, a family of curves C getting closer and closer to the wall; since we need 
estimates uniform in all such curves, one has to be able to deal with such a situation. 
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Figure 13. Left: The curve C; the shaded area represents the convex body whose bound- 
ary is Cn and the dashed line is the boundary of Dj!(iV). Right: The polygonal approxi- 
mation Vn, the dots representing its vertices. 



(in the right order), and staying in some small boxes along its edges. For all edges of 
length smaller than N5n, as well as for the (up to 8) pieces we added above to join 
Cn to the boundary, we impose that the corresponding piece of the contour is a fixed 
length- minimizing path between the vertices. 

The rest of the argument proceeds in a similar way as in the standard case. The 
estimates in the phase of small contours carry over without any problems since in that 
case the effect of the boundary field cannot propagate far away from the wall. 

We still have to explain how one can extract the correct surface tension for Vn from 
the sum over contours in the class (3 introduced above. To do this, we use several results 
about the random-line representation, proved in [ PV2 , PV3 |. To lighten the notation, we 
simply write q N ,v instead of %*''' ; 0* and rf are the dual of /3 and rj, see ( 1.2.5 ). The 

first inequality is just the analogue of ( |3, 4. 6 ) in our case, which turns out to be valid for 
arbitrary ferromagnetic coupling constants: The weight of any high-temperature contour 
7 G © satisfies (|PV|], Lemma 5.4) 



where 7^ denotes the piece of the contour 7 between the fcth and k + 1th vertices of Vn- 
The next step is to replace q N ' (7^) by the corresponding infinite- volume quantity. First, 

for any 7^ joining vertices not belonging to T** N = {i G Ojv.rAT* : *(2) = — \ } (note that 7^ 
stays necessarily at a distance 0(N5n) from T, N ) 

second, for the pieces 7^ joining two sites of H N , we use 



1n' V 



where Lf = {i G Z 



: i(2) ^ — 5} (both results are proved in [PV2], Lemma 5.3). Finally, 

the remaining pieces have a length at most 8N5n, so that their total weight is larger than 

e -CO(NS N )_ 
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The last step is to extract the surface free energy. The basic tool to do this is, as in the 



proof of Theorem 3.4.4, concentration properties for open contours between 2 fixed dual 
sites. For the pieces 7*. not touching the boundary, we can use the usual infinite volume 
results based on ( 5.2.12] ), setting s = N5n- For the pieces along the boundary, one can 



use the following statement ( |pV3 |, Lemma 6.10): 

i""W>W!}'' (1 +«(!)). (4.4.1) 



A: i—>j 
AcN K (i,i)nLf 



where N^(i, j) is defined in Appendix B (with s = N5n)- (In fact, ( 4.4.1 ) can be strength- 
ened when r\ < rj„(P): in this case, the set N^-(i, j)nLf can be replaced by the set ( |P V3 ] , 
Lemma 6.13) 

{fceLf : (i(l)Aj{l))-KlogS N ^k(l) ^ (i(l)Vj(l)) + Klo g 5 N , fc(2) ^ K log 5 N } , 

which is compatible with our picture of partial wetting.) 

The result then follows from lower bounds on the corresponding 2-point functions. The 
only new inputs are the following lower bounds on the boundary 2-point function, 

/ vflw ^ ^ ex p{- r bd(M)IU -»||} w . , m (aai\ 



1 



||3/2 



(Wj)[ d ? > C exp{-r bd (P,v)\\j ~ Vt? < r? w (/3) , (4.4.3) 



for any i,j G X* = {k G Lf : fc(2) = -|}. (gAg ) is proved in gV| Prop. 7.1, while 



(4.4.2) follows from exact computations in the case rj* = 1 | |V1W[ |, and [ PV2| |, Prop. 7.1, 

((Ti^)^"* > (tanh/3*) 2 (a^^f^ 1 , Vr? ^ . 



Second case: 77 = 

This is a somewhat marginal case. The apparent difficulty is that in this case rf = 
00. However, this does not create any real complications. One just has to modify the 
construction of the first case as follows: We replace the polygonal line Vn by the (possibly 
open) polygonal line Pn \ {u G M 2 : n(2) =0}; we then sum over contours going through 
the vertices of this polygonal line (contours which are open if the polygonal line is open). 
This does not give any contribution for the part of C along the wall, which is what we 
want since r bd (/3, 0) = 0. 

Third case: n < 

This is slightly more tricky. In this situation, one may be even more pessimistic, since the 
duality is simply not defined when non- ferromagnetic interactions are present! However, 
this turns out to be a false problem. Indeed, we can use the following obvious identity to 
recover ferromagnetic interactions (see footnote [j| p. 53), 



P,ri _ AM 

where ± correspond to the boundary condition a, = 1 if i{2) ^ and <fj = — 1 otherwise. 

We then construct Vn as in the first step and set / = Vn H {x G B 2 . : x(2) = 0}. 
If / = 0, then we subdivide the set {x G O 2 : x(2) = 0} into segments of length <5jv/2 
(except possibly for the last one, which might be shorter); this defines a second (open) 



polygonal line V' N (with all its vertices along the wall) (see Fig. 14). We then introduce a 
class of pair of contours (7, 7'), 7 going through the vertices of Vn and defined as before, 
and 7' following the wall, going through the vertices of V' N and staying inside small boxes 
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Figure 14. The construction for r\ < 0. Left: 1 — (two polygonal lines: one open and 
one closed. Right: 1^0 (one open polygonal line). 



along its edges, similarly as for the other one (7' is open). By construction 7 and 7' are 
disjoint. Duality then implies the following identity 

M5g ({7, 7'} C 7 ( • )) = (Z©- 1 ffrMV) E 





(Ci7i7') A*-comp. 



= (l-- 0(W) )^4' H (7,y). (4.4.4) 

Z N,± 

The factor (1 — e~ oi - N ^) comes from the fact that we can apply duality only to simply 
connected sets, and the exterior of 7 is not simply connected. We must therefore forbid 
families £ for which duality does not hold; since such families must contain at least one 
contour surrounding 7, we get the above correction. 

We can now proceed as in the first case. The only additional work to do is to analyze 
the ratio of partition functions in ( [4.4.4] ), but this is easy, since by duality 

N >+ - (l a . a . > r hd (J3,\ V \)(2N+l) u 4 ^ 

where t\ = (— L — 5, — |) and t x = (L + — h) are the two dual sites at the lower left and 
lower right corners of W N rN , and the last inequality follows from the upper bound (see 



|PV2| for example) 



valid for any i,j £ S^-. We then see that the ratio of partition function cancels the 
contribution from the sum over the open contour 7', up to an error term exp{0(N5]\f)}. 

If I ^ 0, the situation is simpler. Let's write I = [x, y}; then we define a new polygonal 
line V^: goes from the lower right corner of to a along the wall, then it follows 
"Pn \ {1 £ D, : x{2) = 0} up to b and finally goes from b to the lower right corner of 
(see Fig. [L4|). We subdivide as usual the part of along the wall into segments of length 
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Figure 15. Left: A contour touching the wall and the family (t>i, . . . , V2m)- Right: An 
s-skeleton for the contour. 



5 n/2 and proceed as in the first case, with replacing Vn-, using ( 4.4.4j ). Summing over 



the open contour going through the vertices of produces (up to the usual error term) 



a term exp{— W, 

exp{2r bd (/3, \ V \) N} eKp{-W m {P%) N} 
the conclusion follows as in the usual situation. 



AtolC^iv) -^J" - Combining this with ( 4.4.5| ) and observing that 



exp{-W/3, v (V N ) N} 



The upper bound. Let us now turn our attention to the proof of the upper bound. The 



basic strategy is completely similar to that of the standard case, see Subsection |3.5.2| . The 
only serious modification concerns the energy estimate, which should now associate the 
functional W/3 )7? to the probability of skeletons. Again, the case rj ^ is somewhat simpler 
than the other, so we start with this one. 

First case: ij ^ 

The basic problem we encounter when trying to make the energy estimate is the same 
we met in the proof of the lower bound. Summing over an open contour connecting two 
dual sites i and j might not yield a decay of order exp{—Tp(J—i)} or exp{— r bd (/3, 77) || j — i\\} 
if i and j are close enough to the wall but not on it (see |PV3| ] ) . However, the following 
bound, proven in | PV2fl , Lemma 5.1, is sufficient to derive the energy estimate, 



A: i— fj 



In 



(A) < exp{-r/3(j - i)} , 



(4.4.7) 



for any r? 0; 6 (Str) = {e* C S^}. The definition of skeletons will be done in such a way 
as to ensure that the additional constraint A Pi £ (S^) = is automatically satisfied, see 
below. We also need to extract the wall free energy when summing over contours joining 
two dual sites belonging to S^; this however is nothing else as ( 4.4.6| ). 

Let us now describe the construction of a skeleton S = (ui, . . . , u n ) of a closed contour 
7. Remember that we have to define the skeletons in such a way as to ensure that 1) the 
piece of the contour between two dual sites not both on the wall must be edge-disjoint 
from the wall, and 2) the Hausdorff distance between the contour 7 and the polygonal line 
Pol(S') is smaller than the cutoff parameter s(N). 

For contours 7 which do not touch the wall, the definition of skeletons is the same as 
in Part ||. Suppose 7 n / 0- Let us define (v\, . . . , V2 m ) as the minimal family of 

dual sites satisfying the following properties: 
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Figure 16. Left: The family of polygonal lines associated to S . Right: The family of 
closed polygonal lines associated to 6. 

1. Vk G Yj* n Pi 7 for k = 1, . . . , 2m and Vk(l) < ffc'(l) if k < k 1 ; 

2. (vi,..., v m ) split 7 into pieces 71 : v\ — > v 2 , ■ ■ ■ , 72m : v 2m — > v\, such that 

• 72fe H <?(S^) = for all k = 1, . . . , m. 

• %(72fc,{^ e E 2 : x(2) = -1/2}) > s(N) for all A; = 1, . . . , m. 

• dmll2k+i,{ x e E 2 : x(2) = -1/2}) < s(JV) for all A; = 1, . . . , m. 
We then say that S = (u±, . . . , u n ) is an s-skeleton of 7 if 

• All vertices of S belong to 7. 

• Vi, . . . , V2m are vertices of S. 

• The only vertices of S along 72^+1 are V2k+i and v 2 k+2, for all k = 1, . . . , m. 

• The distance between any successive pair of vertices ui,ui + \ of 5 along 72ft satisfies 
s(N)/2 < Ih-uj+ilU < 2s(iV), for all A; = 1, . . . , m. 

• d H (7,Pol(5)) ^s(N). 

This definition has the nice property that either ui and ui+\ both belong to S^, or the 
part of 7 between these two sites is edge-disjoint from T,* N (see Fig. 15). This allows us to 



use the estimates (|4.4.6|) and Qi.4.7| ) . This yields the following extension of fl3.4.2|) JPV2 

4?* (©)< exp{-W^(6)}. (4.4. 



The analogue of the energy estimate ( 3.4.3| ) then follows easily, since r bd (/3, rj) ^ when 



rj ^ and therefore it is still possible to control the number of vertices of S in terms of 
W/3 )J7 (6). This gives 

Using this and the estimates in the phase of small contours, which still hold in the presence 
of a boundary field, the upper bound follows easily. 

Second case: rj < 

As for the lower bound, we have to deal with the fact that, for rj < 0, the duality is not 
defined. The solution is the same as there: We just change boundary conditions, i.e. we 
look at the measure /ijy + , which was defined when we dealt with the lower bound. 

Once we have done this, the main difference is that the family of low-temperature con- 
tours of any configurations compatible with these boundary conditions contains exactly 
one open contour, with endpoints t\ = (-N — \, — |) and t r = (N + ^, ^). It is straightfor- 
ward to generalize the notion of skeleton introduced in the preceding case to the present 
situation. What we get by this procedure is a family of skeletons = (So, Si, ... , S n ) 
containing exactly one skeleton, So, with PoI(Sq) open with endpoints t\ and t r . 
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Since we want to compare the corresponding families of polygonal lines with the solution 
of the variational problem, i.e. with the boundary of a convex body in D;5, it is convenient 
to introduce another family 6 of skeletons whose associated polygonal lines are closed; & 
possesses the same set of vertices (except for t\ and t T , but with a different set of edges, 
which is such that its associated family of polygonal lines satisfies 

Pol(S) = Pol(S ± )A{x G M 2 : -N/2 - \ < x(l) < N/2 + ±, x{2) = -±} 

where A denotes symmetric difference (see Fig. [l6| ). 
One then has the following relation 

W A „(6) = W AW (e±) - (2N + 1) r bd (/3, \r)\) . 

In particular, the following version of ( (4.4.8|) holds |PV2| 1 

/ijg(6 ± ) ^ if! exp{-W^(S)} r, > - V M (4.4.10) 

/4±( S± ) < K 2 N 3 / 2 exp{-W^(S)} r, < - r? w (/3) (4.4.11) 



The energy estimate ( [4.4.9 ) is slightly more delicate now, since the wall free energy is 



negative. It turns out however that in the partial wetting regime, rj > — r/ w (/3), it is easy 
to reduce ourselves to a situation similar to the case r bd (/3, rj) ^ 0. The case rj ^ — r/ w (/3), 
i.e. complete wetting, is more subtle, but happens not to give too much problems as 
long as we consider volume-order large deviations (or, in fact, deviations close enough to 
volume order). 

Let us first consider the case of partial wetting; this regime is characterized by |r bd (/3, r/)\ < 
Tg. Let us write W/? )r) (<5) = T + + T~, where T + (T~) is the positive (negative) part of 
the functional. Then, since T + ^ {r% / T hd (fl , rj)) T~ and the number of vertices along the 
wall is at most two-third of the total number #(6), we have 

for some absolute constant K. This allows to prove that 

/$L(6 > r) < exp{-r(l - ^^)} • (4.4.12) 

When rj ^ — r/ v (f3), one cannot establish so good an upper bound. The best we can do is 
to use the fact that T~ ^ (2N + 1) r bd (/3, rj), which turns out to be enough to prove the 
following, weaker, version of the energy estimate 

/^(6 r) < exp {- r( l - + C M} . (44.13) 

The reason why such an estimate is still sufficient to get the desired result is that the 
relevant values of r are also of order N, so that the first term can always be made to 
dominate the second one. 



Once we have (4.4.12) and (4.4.13), the proof is concluded as usual, after observing that 



the estimate in the phase of small contours still applies in the presence of the boundary 
field \rj\. 
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4.4.2. Ising model in D ^ 3. The proof of Theorem 4.3.3 is based on the Li-Theory 
introduced in Part [|. We simply explain how the main ingredients of the proof should be 
modified and refer to |BIV] for details. 



The arguments of geometric measure Theory can be extended easily to this new set- 
ting. In particular, it is straightforward to check that the functional W/3 j?? is lower semi- 
continuous and that the approximation Theorems p. 5.1 and 2.6.1 hold. 

The main problem is to define proper mesoscopic phase labels for the measures with a 
boundary magnetic field. If rj ^ 0, then the mesoscopic phase labels introduced in Part ^ 
satisfy the Assumptions A and B, as well as Conditions C1-C3 under the measure 
Instead if rj < 0, some problems occur because the FK measure looses its ferromagnetic 
properties and the random coloring measures are more complicated to deal with. Never- 
theless, it is still possible to define mesoscopic phase labels and to derive estimates as in 



Section 2.2. 



Other difficulties have to be overcomed in order to implement the general philosophy of 
the Li-Theory. In the case of a negative boundary magnetic field, the interface induced 
by the field prevents us from applying directly the techniques developed to prove the ex- 
ponential tightness Theorem 2.1.1| . Therefore an alternative approach similar to the one 
described in Subsection [4.4.1 is required. The analysis of the surface tension needs also 



some care. We recall that the computation of surface tension is based on a localization 
procedure along the boundary of functions of bounded variation. For a given test function 
either locally its boundary is in the bulk and we recover the usual surface tension term or 
it intersects the wall and arguments similar to those used in the bulk enable us to derive 
the wall free energy. In this way the complexity of the problem is reduced because the 
difficult analysis of the fluctuations of the microscopic interface between the wall and the 
bulk is replaced by soft Li estimates. 



4.5. Open problems 

As in the previous parts, there are still a lot of open problems. Most of those pre- 
sented before have natural analogues in the present situation. In the following, we restrict 
ourselves to problems intrinsically related to the topics discussed in this part. 

2D nearest-neighbors Ising model. The fact that one is still unable to analyze non-pertur- 
batively the fluctuations of the phase separation line is only strengthened when we would 
like to study boundary effects. Indeed, a general analysis of typical open paths with 
endpoints at general positions with respect to the wall has not been done even at low 
temperature. Problems related to this are the following: 

1. Give a non-perturbative proof that the probability measure of a suitably rescaled ver- 
sion of an open contour with endpoints on the wall converges weakly to the measure 
of Brownian excursion when r\ ^ — rj w (/3) (as was sketched in the low-temperature 
case for rj = —1 in 0). This would provide a way of analyzing the typical fluc- 
tuations of magnetization in the complete wetting regime, and would complete the 
heuristic picture of the wetting transition in the Grand-Canonical Ensemble. 

2. Establish Ornstein-Zernike behavior for the boundary 2-point function without hav- 
ing recourse to explicit computations. Even weaker lower bounds, like those given in 
[A] |, have not been proved in such a constrained geometry. 



OS") 



Another open problem is to investigate the full range of moderate deviations. This may 
require an understanding of point 1. above. 

Higher dimensional nearest-neighbors Ising models. If fluctuations of phase separation 
lines are not yet understood, the situation is only much worse when considering their higher 
dimensional counterparts; in fact, even perturbative results are not always available. Here 
is a far from exhaustive list of related open problems. 

1. Give a microscopic description of the behavior of phase boundaries in the partial and 
complete wetting regimes in the Grand- Canonical Ensemble to put some flesh on the 
heuristics given above. 

2. Decide whether t? w (/3) = 1 or not. The corresponding results for the SOS model Cb | 
suggest that 7/ w (/3) < 1 in any dimension; numerical investigations confirm this in 
dimension 3 |BL], 

In fact, even much simpler problems related to behavior of higher dimensional interfaces 
are still open: proof of the existence of a roughening transition in d = 3, proof of the 
unstability of the (1, 1, 1) interface, ... 

In some simpler models of the SOS type some (but not all!) of these problems can be 
solved, but this does not seem to help in solving the original ones. 

The wall. Another type of problems concerns properties of the wall. In particular, it might 
be interesting to answer the following questions. 

1. What happens if the interaction with the wall is more complicated (say, non-nearest 
neighbor). 

2. What happens if the boundary field is not homogeneous (for example, is a "random" 
configuration of r\\ and 772 macroscopically equivalent to some well-chosen homoge- 
neous boundary field r] = ff?). 
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Part 5. Appendix 



5.1. Appendix A : Proof of Theorem 2.2.1 



Assumption A controls the number of zero u^-blocks, whereas Assumption B is used to 
control the geometry of the mesoscopic phase labels. The dependence of ko on 5 could be 
described as follows: we choose ko so large that 

1 

Pk < 7773T<5 for ever y k ^ (5.1.1) 

C(d) 

where C(d) is a large enough fixed constant. Three terms on the left hand side on ( |2.2.3| ) 
correspond to three different exponential estimates: 

5.1.1. Estimate on the volume of zero u^-blocks. The domination by Bernoulli mea- 
sure ( [2.2.2 ) implies that 



N 



N\ d \ { /N\ d . 5 



Each realization of the phase label function Uk splits T d into the disjoint union of three 
mesoscopic regions: 

f d = {x: u k (x) = 1} V {x : u k (x) = -1} V {x : u k (x) =0} = A + V A_ V A . 

By the choice of the scale ko in (|5.1.1[) the estimate ( |5.1.2| ) is non-trivial for every k ^ ko, 
and, in view of the target claim fl2.2.3j ) , we can restrict attention only to such realizations 
of u k for which 



|A | = / l{ Uk ( x )=o}&x < 5. (5.1.3) 

This has the following important implication: if u k 6 V (K a , 25) c , the area of the boundary 
of any regular set A such that A + C A C T d \ A_ is bounded below as 

\dA\ ^ a. (5.1.4) 

Using the Assumption B of the Theorem we are going to construct such sets A on the 
finite ko scale; A £ ^F n -ko > an d m sucn a fashion that all the boundary fco-blocks of A will 
necessarily have zero u^-labels. This reduction enables a uniform treatment of all coarser 
scales k ^ ko. 

So let k ^ ko, and assume that ( |5.1.3 ) holds. We denote by j4_ (resp. A + ) the set of all 



boxes B n _fc in A_ (resp A + ). We say that x £ T^_ fco is — * connected to j4_; x < — -> A_, 

if there exists a *-connected chain of "— " u ko blocks leading from M n _ ko {x) (and including 
it) to A^. Define now the complement A c as follows: 

A c = A^ |J B n _ fe0 (x). 

By the virtue of the Assumption B, A + C A. Moreover, by construction all the fco-blocks 
of A attached to the boundary dA° have zero Ufc -labels. With a slight abuse of notation 
we proceed to denote this collection of boundary /co-blocks as dA. By ( 5.1.4| ) the number 
of fco-blocks in dA is bounded below by 

^o(OA) >^N*~\ (5.1.5) 
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Since, however, the total number of fco-blocks i n the corresponding decomposition of 
equals to N d /2 dk ° the estimate fl5.1.5| ) alone is not sufficient for giving the desirable upper 
bound on the probability P^v (uk £ V(iT a , 2<5) c ). The required entropy cancelation stems 
from the fact that small connected contours of dA cannot surround too much volume. 
Let us decompose A to the disjoint union of its maximal connected components: 

i i 
A = \l Ai respectively dA = \l dAi. 
i=\ i=i 

We shall quantify contours dAi according to the size (or the number of feo-blocks ) in Ai. 
Namely, the contour dAi is called small, if 

# ko (Ai) ^K(d)logN or 1^1 ^K(d)^logN, (5.1.6) 

where K(d) is a sufficiently large constant. Otherwise, the contour dAi is called large. 
We claim that under ( 5.1. 3D the following inclusion is valid: 



{u k GV(K a ,25) c } C { 



E MMx^Ui E \ dA i\> a \- (5-1-7) 

dAi— small I I dAi — largo 



Indeed, if the total volume inside small contours is less than 5, then repainting all the 
small components Ai into "—1" and all the large components Aj into "+1" we produce a 
{±l}-valued function which is at most at the Li-distance 25 from u k and which, thereby, 
cannot belong to K a . 



5.1.2. Peierls estimate on the size of large contours. 



p* I E \°m >*)=**[ E mm) > ^k^- 1 , 

\9yli-large / \dAi-large J (5.1.8) 

This immediately follows from Assumption A, once the constant K(d) in ( 5.1.6| ) has been 
properly chosen. 



5.1.3. Estimate in the phase of small contours. The volume of small components Ai 
is related to the total number of fco-blocks in these components as 

-d 

dAi— small dAi— small 



E w = (£) " E #*o(Ai 



On the other hand, for every I G [1, n — ko]; 

E #k (Ai) = ^2 E 1-^eA;} 

dAi— small x£T d dAi— small 



E E E l{f*A a;6Ai} ; 

te[0,...,2') d xe j ^„ fco _ ( 9^i-small 
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where Ao = 2 k ° n is the step size on the embedded torus T^_ fco , and 9, is the shift on 
this torus. Consequently, 

E W > 5 ) < te K/"( E E W^}>4^n)' 

\dAi— small J yaiGT 1 * fc ; 9Ai — small 

(5.1.9) 



If, however, 2 l > if (d) log N, then no two distinct points on the torus T^_^_j (or any 
shift of it) can belong to the same small component A{. This, in view of the domination 
by the independent Bernoulli site percolation (Assumption A), suggests an application of 
the B-K inequality. Since, by the choice of the scale &o i n ( |5.1.lD ; 

£fc = Pperc (3 a closed surface of zero li/^-blocks around x) < 5, 



for every x E , we readily obtain that the right hand side of ( 5.1 .9| ) is bounded above 

by 



c 4 (d)exp { -5 ( ) log ( — 



e k . 



The proof of Theorem 2.2.1 is concluded. □ 



5.2. Appendix B : Proof of the three-point lower bound Lemma 3.4.3 



The proof of Lemma 3.4.3| is based on the following positive stiffness property of the 



surface tension [AA|: 



{ Jr" ( " (9)> + T " im) } = 4« " 9 im) > a <5 ' 2 ' 10) 

where the unit normal n{9) is defined via n{9) = (cos 9, sin 9), and Rp (n) is the radius 
of curvature of dK, at the point supporting the tangent line orthogonal to n. An integral 
version of ( |5.2.10| ) is the strong triangle inequality |V| : For any «,ti£l 2 : 

T/3 (u) + T0 (v) - T/3 (u + v) > ci(/3) (||u|| 2 + \\v\\ 2 - \\u + u|| 2 ) . (5.2.11) 

The latter inequality is used to control the fluctuations of the microscopic phase boundaries 
(in their random line representation of Section |3.4[) . 

Let now an (s, e)-compatible triple of points (u, u>, v) be given. Fix K = K{[3) large 
enough and define the "oval" neighborhood N^(n, w) of {u, v} as: 

Nk{u,w) = {z E M 2 : Tp (z — u) + T/3 (w — z) — Tp (w — u) ^ K logs}. 

The oval neighborhood N;{(tt),)j) is defined exactly in the same fashion. Relations fl3.4.8| ) 
and (|3,4.11| ) readily imply that that the main contribution to (<J u &w) r (respectively to 



(a w a v )^ ) comes from the paths Ai (respectively A 2 ) which stay in N^(ti, w) (respectively 
Nk(w,v)). More precisely, 



Q P * (Ai) > K^)f (1 + o(l)) , (5.2.12) 



AeN x (u,m) 
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uniformly in all (s, e) -compatible triples. Any such path Ai = (Ai(0), Ai(ni)) could be 
decomposed as follows: Define 

n w = max {A; : A& G Nk(u,w) \ Nk(w,v)} , 

and set A" = (Ai(0), Xi(n w )), Xf = (Xx(n w + 1), Ai(m)); Ai = A" V X'f. The de- 
composition A2 = Af V A?f is defined in a completely symmetric way. Notice that, by the 
construction, the paths A" and A?> are disjoint and compatible, and, by ( 5.2.11] ); 

max{||Ai(n u ') - w\\ 2 , \\X 2 (n w ) - io|| 2 } < c 2 (e) logs. 

The claim of the lemma follows now from Q3.4.6 ) and ( |3.4.7| ), □ 
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